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Abstract 

Aquatic  cohesive  sediments  are  often  the  final  receptor  of  contaminants 
released  by  Department  of  Defense  (DoD)  activities.  In  estuaries,  harbors, 
and  channels,  fine  sediments  can  be  mobilized  and  transported  by  the 
combined  action  of  both  relatively  steady  currents  (forced  by  river  flow, 
tides,  and  wind)  and  by  the  larger  unsteady  bottom  shear  stresses 
associated  with  higher-frequency  oscillatory  wave  motions  (either  wind-  or 
vessel-generated  waves).  Methods  to  measure  contaminated  sediment 
stability  and  erosion  as  influenced  by  wave  processes  are  critical  to 
assessing  contaminant  fate  and  associated  risk  as  well  as  evaluating 
remediation  options.  This  research  focused  on  calibrating  and  verifying 
the  applicability  of  a  prototype  flume  designed  to  measure  erosion  of  fine¬ 
grained  cohesive  sediment  under  combined  unidirectional  and 
superimposed  oscillatory  bottom  shear  stress.  These  are  similar  to  the 
shear-stress  conditions  that  exist  in  wave-current  environments.  This 
research  will  provide  more  accurate  methods  for  assessing  contaminated 
sediment  stability  for  many  DoD  and  Environmental  Protection  Agency 
managed  contaminated  sediment  sites. 
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Executive  Summary 

Objective 

The  proposed  research  goal  is  to  develop  laboratory  methods  to  reliably 
measure  and  predict  erosion  of  contaminated  fine-grained  sediments 
under  bottom  shear-stress  conditions  that  represent  combined  wave  and 
current  action.  The  shear-stress  conditions  for  erosion  measurements  will 
replicate  those  expected  in  the  field  (nonbreaking  wave,  shallow  water, 
estuarine,  harbor,  and  coastal  environments).  The  specific  objective  of  this 
project  is  to  calibrate  and  validate  an  innovative  mobile  erosion  flume 
designed  to  measure  cohesive  sediment  erosion  under  combined  wave- 
current  conditions. 

Background 

Aquatic  cohesive  sediments  are  often  the  final  receptor  of  contaminants 
released  by  Department  of  Defense  (DoD)  activities.  In  estuaries,  harbors, 
and  channels,  fine  sediments  can  be  mobilized  and  transported  by  the 
combined  action  of  both  relatively  steady  currents  (forced  by  river  flow, 
tides,  and  wind)  and  to  an  even  greater  extent  by  the  larger,  unsteady 
bottom  shear  stresses  associated  with  higher-frequency  oscillatory  wave 
motions  (either  wind-  or  vessel-generated  waves).  Methods  to  measure 
contaminated  sediment  stability  and  erosion  as  influenced  by  wave 
processes  are  critical  to  assessing  contaminant  fate  and  associated  risk  as 
well  as  evaluating  remediation  options. 

Materials  and  methods 

This  research  focused  on  calibrating  and  verifying  the  applicability  of  a 
prototype  flume  designed  to  measure  erosion  of  fine-grained  cohesive 
sediment  under  combined  unidirectional  and  superimposed  oscillatory 
bottom  shear  stress.  These  are  similar  to  the  shear-stress  conditions  that 
exist  in  wave-current  environments.  The  new  flume  is  based  on  the  design 
of  a  well-documented,  frequently  utilized  flume  for  measuring  cohesive 
sediment  erosion  under  unidirectional  bottom  shear  stress  representing 
steady  current  (Roberts  et  al.  1998;  McNeil  et  al.  1996).  Preliminary  results 
using  the  prototype  flume  (Sediment  Erosion  Actuated  by  Wave 
Oscillations  and  Linear  Flow  [SEAWOLF])  for  the  combined  unidirectional 
and  oscillatory  flow  have  been  published  (Jepsen  et  al.  2002). 
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Bottom  shear  stress  under  wave/ current  conditions  has  been  studied 
extensively  in  the  field  and  in  large-scale  laboratory  wave  flumes  (Grant  and 
Madsen  1979;  Dibajnia  and  Watanabe  1992;  Wikramanayake  and  Madsen 
1994).  The  time  history  of  shear  stress  that  the  sediment  bed  experiences 
under  these  combined  conditions  is  therefore  well  understood  in  theory, 
and  algorithms  exist  that  accurately  predict  shear  stress  (Bijker  1971;  Swart 
1976).  The  SEAWOLF  flume  does  not  replicate  the  wave/current  hydro- 
dynamic  conditions  but  rather  replicates  only  the  bottom  shear  stress  that 
causes  erosion.  Free  surface  waves  creating  high-bottom  shear  cannot  be 
developed  in  a  field-deployable  device  at  reasonable  cost.  The  interior  of  the 
flume  channel  must  be  instrumented  to  demonstrate  that  bottom  shear 
stress  over  the  sediment  surface  in  the  flume  is  equivalent  to  specified  wave- 
current  conditions.  The  first  step  in  demonstrating  flume  validity  to  wave- 
current  conditions  is  to  demonstrate  that  bottom  shear  stress  over  the 
sediment  surface  in  the  SEAWOLF  flume  replicates  shear  stress  due  to 
combined  wave-currents.  To  demonstrate  this,  two  instruments  have  been 
incorporated  into  the  flume,  one  for  measuring  both  near-bed  and  full-field 
flow  structure  (including  velocity  profile)  and  the  second  that  simply 
measures  near-bed  velocity  profile  from  which  shear  stress  can  be  derived. 
The  first  system  is  a  novel  high- magnification,  particle-based  velocimetry 
system  that  utilizes  particle-image  velocimetry  (PIV)  and  particle-tracking 
velocimetry  (PTV)  algorithms  to  resolve  micro-  and  macro-scale  flow 
structure  as  well  as  the  velocity  profile  to  within  7  /an  of  the  test  section 
wall.  Classical  velocimetry  systems  are  typically  deployed  in  lower-energy, 
steady-flow,  free-surface  environments.  Implementation  within  a  closed- 
channel,  high-shear,  rapidly  time-varying  environment  is  a  new  application 
for  these  devices  and  one  of  the  most  challenging  research  aspects  of  this 
project.  Issues  related  to  time  scale,  turbulence,  and  secondary  circulation 
under  which  these  instruments  have  not  been  tested  had  to  be  analyzed. 

The  second  system,  the  miniLDV,  is  a  self-contained  and  permanently 
aligned  laser  Doppler  velocimeter.  The  miniLDV  probe  contains  a  laser, 
miniature  beam  shaping  optics,  receiving  optics,  and  a  detection  system. 
The  size  of  the  probe  volume  (measurement  domain)  is  -25  /an  by  50  /an 
by  150  /an.  The  system  is  mated  with  a  computer-controlled  linear 
traversing  system  allowing  for  resolution  of  near-bed  velocity  profiles. 


Based  upon  the  very  interesting  results  from  the  hydrodynamic  data 
collection  activities,  and  with  the  Strategic  Environmental  Research  and 
Development  Program  (SERDP)  review  panel  approval,  the  path  forward 
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changed  in  2008  to  focus  exclusively  on  hydrodynamic  calibration  and 
validation.  Verification  is  required  to  prove  that  the  bottom  shear  stresses 
created  in  the  flume  represent  bottom  shear  stresses  under  naturally 
occurring  wave-current  conditions.  Demonstrating  flume  applicability  to 
wave-current  environments  requires  designing  standard  operating 
procedures  for  SEAWOLF  flume  operation  and  collecting  sufficient 
comparative  laboratory  data  to  prove  SEAWOLF  shear-stress  replication 
methodology. 

The  addition  of  the  miniLDV  sensor  will  bring  significant  added  benefit  to 
the  project,  providing  a  coordination/comparison  of  several  state-of-the- 
art  instruments  never  before  tested  in  this  type  of  environment.  This 
rapidly  varying,  enclosed,  high-shear  environment  is  required  not  only  for 
this  project  but  for  all  cost-effective,  high-energy  sediment  erosion 
experiment  methods.  Cost-effective  methods  will  require  replicating  shear 
stress  instead  of  hydrodynamic  conditions.  Therefore,  the  benefits  of  the 
shear-stress  validation  and  procedure  comparison  will  greatly  benefit  the 
larger  sediment  transport  research  community. 

Results 

Two  major  hydrodynamic  studies  are  presented  in  this  document.  First, 
PIV  was  used  to  study  the  cyclic  modulation  of  the  wall  shear-stress  and 
turbulence  properties  of  an  oscillating  channel  flow.  The  PIV  instrument 
employed  utilized  a  dynamically  adjusted  delay  between  the  laser  pulses  to 
accommodate  the  wide  variations  in  velocity  encountered  in  the  oscillating 
flow.  Both  high-  and  low-magnification  digital  PIV  recordings  were 
obtained  to  reveal  the  near-wall  boundary  layer  structure  and  wall  shear 
stress,  as  well  as  the  full-field  turbulence  throughout  the  channel. 
Presented  are  wall  shear-stress  and  global  turbulence  data  for  Stokes- 
thickness  Reynolds  numbers  of  Res  =  1220,  2033,  and  2875.  The  results 
reveal  a  fully  developed  turbulent  state,  relaminarization,  and  an  explosive 
transition  back  to  turbulence.  The  flow  is  examined  in  detail  for  the  case  at 
Res  =  1220,  where  instantaneous  PIV  realizations  at  low  magnification 
reveal  the  structure  of  the  flow  during  relaminarization  and  transition 
back  to  turbulence.  High-magnification  PIV  results  are  used  to  reveal  the 
phase  modulation  of  the  mean  velocity  profiles  in  the  viscous  sublayer  and 
logarithmic  layers  through  the  half  cycle  and  quantitative  profiles  of  in¬ 
plane  Reynolds  stresses  and  turbulence  production  are  presented.  To  our 
knowledge,  this  is  the  first  PIV  investigation  of  this  canonical  unsteady 
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turbulent  channel  flow,  and  these  results  represent  a  needed  contribution 
to  the  limited  turbulence  data  that  exists  for  unsteady  wall  flows. 

Second,  the  structure  of  turbulence  in  an  oscillating  channel  flow  with  near- 
sinusoidal  fluctuations  in  bulk  velocity  is  investigated.  Phase-locked 
particle-image  velocimetry  data  in  the  streamwise/wall-normal  plane  are 
interrogated  to  reveal  the  phase-modulation  of  two-point  velocity 
correlation  functions  and  of  linear  stochastic  estimates  of  the  velocity 
fluctuation  field  given  the  presence  of  a  vortex  in  the  logarithmic  region  of 
the  boundary  layer.  The  results  reveal  the  periodic  modulation  of  turbulence 
structure  between  large-scale  residual  disturbances,  relaminarization 
during  periods  of  strong  acceleration,  and  a  quasi-steady  flow  with  evidence 
of  hairpin  vortices  that  is  established  late  in  the  acceleration  phase  and 
persists  through  much  of  the  deceleration  period.  Again,  to  the  best  of  the 
authors'  knowledge,  this  is  the  first  structural  information  for  an  oscillating 
flow  derived  from  direct  velocity  measurements. 

Also  in  2008,  a  newly  developed  commercial  miniLDV  sensor  was  identified 
for  potential  application  within  the  SEAWOLF.  While  reliable,  PIV-derived 
measurements  of  shear-stress  time  histories  offer  limited  amounts  of  data. 
Characterization  of  a  single  wave-current  frequency  and  amplitude 
operating  condition  generally  requires  a  full  day  of  testing  and  an  additional 
1-2  days  for  analysis  of  the  PIV  data  for  the  extraction  of  wall  shear-stress 
histories.  For  the  purposes  of  obtaining  large  parametric  data  sets  that  vary 
in  frequency,  amplitude,  and  bed-roughness  conditions,  a  more  rapid  shear- 
stress  measurement  must  be  incorporated.  In  previous  years,  a  floating 
element  mechanical  shear  stress  sensor  was  tested  that  could  be  calibrated 
to  a  shear-stress  value.  These  mechanical  sensors  have  been  applied  for 
testing  in  SEAWOLF  with  mixed  results;  the  sensors  responded  to  time- 
varying  shear  stresses  but  displayed  weak  signal  levels  that  were  easily 
corrupted  by  electrical  noise  and  suffered  from  time-dependent  drift  in 
sensor  voltage,  which  corrupted  the  measurement  as  well. 

As  an  alternative  to  mechanical  shear-stress  measurements,  two  newly 
developed  commercial  optical  sensor  technologies  were  tested,  the  microS 
and  miniLDV,  which  (1)  offer  calibration-free  measurements,  (2) 
improved  signal-to-noise  characteristics,  and  (3)  have  been  proven  to 
work  in  harsh  and  dirty  environments.  The  miniLDV  significantly  out¬ 
performed  the  microS  in  initial  testing,  and  therefore  it  will  be 
implemented  within  the  SEAWOLF  facility. 
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Benefit 

The  proposed  research  will  provide  more  accurate  methods  for  assessing 
contaminated  sediment  stability  for  many  DoD  and  Environmental 
Protection  Agency  managed  contaminated  sediment  sites.  The 
hydrodynamic  validation  exercises  will  produce  stand-alone  products  that 
will  benefit  multiple  other  sediment  transport  projects.  The  final  methods 
developed  for  measuring  and  predicting  contaminated  sediment  erosion 
under  wave-current  conditions  will  be  peer  reviewed,  cost  effective,  and 
standardized.  Such  methods  do  not  presently  exist,  adding  significant 
uncertainty  in  predicting  contaminant  fate  and  associated  risk  in  shallow- 
water,  coastal,  estuarine,  and  harbor  environments.  The  proposed  research 
will  develop  methods  to  reduce  uncertainty  that  presently  complicates  risk 
assessment  and  evaluation  of  remediation  options. 

Transition  plan 

There  is  a  strong  need  for  capabilities  to  measure  erosion  rates  for  buried 
contaminated  sediment  under  wave-current  conditions.  A  similar  need 
previously  existed  for  steady  (current  only)  conditions.  The  transition  of 
SEDflume  (the  current-only  equivalent  to  the  proposed  wave-current 
flume)  from  research  tool  to  frequently  used  field  device  for  contaminated 
site  evaluation  was  rapid  (fewer  than  3  years)  due  to  the  strong  field-user 
need  and  well-documented,  peer-reviewed  flume  validation.  Similarly,  the 
transition  of  the  SEAWOLF  flume  from  research  tool  to  site  evaluation 
device  is  expected  to  be  rapid  for  two  reasons:  (1)  the  flume  will  be 
validated,  and  there  are  no  alternative,  cost-effective  field  devices  for 
measuring  erosion  rates  under  these  wave-current  conditions  and  (2) 
there  is  a  strong  need  in  the  remediation  community  for  quantifying  wave- 
current  erosion  of  fine-grained  sediments.  An  additional  transition  phase 
will  involve  incorporating  erosion  algorithms  into  existing  predictive 
models  such  as  the  Advanced  Hydraulics  Model  (AdH)  and  Long  Term 
FATE  (LTFATE)  model  (U.S.  Army  Engineer  Waterways  Experiment 
Station  1998). 
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1  Objective 

This  Strategic  Environmental  Research  and  Development  Program 
(SERDP)  Project  addresses  SERDP  Statement  of  Need  CUSON-06-03, 
Assessment  and  Measurement  of  Processes  Impacting  the  Fate  and 
Transport  of  Contaminants  in  Sediments.  The  project  specifically  addresses 
the  following  high-priority  needs  outlined  in  the  Statement  of  Need  (SON): 

Develop  and  evaluate  site  characterization  tools  to  measure  the 
rates  of  important  sediment  physical  processes  affecting  the  fate 
and  transport  of  contaminants,  and  develop  improved  methods  for 
incorporating  uncertainty  into  measurements  of  fundamental  fate 
and  transport  processes. 

The  overarching  goal  of  Project  ER 1497  is  to  develop  methods  to  reliably 
measure  and  calculate/predict  erosion  of  buried,  contaminated,  fine¬ 
grained  sediments  under  combined  wave  and  current  stresses.  The 
methods  will  be  developed  for  wave-current,  shear-stress  regimes  typically 
encountered  at  contaminated  sediment  sites  (shallow  water,  harbor, 
estuarine,  and  coastal  environments).  The  measurement  devices  will  be 
housed  in  a  mobile  trailer  so  that  the  measurements  can  be  obtained  at  the 
site  with  minimum  disturbance  to  the  tested  sediments. 

This  project  includes  the  following  technical  objective: 

Calibrate  and  validate  an  innovative  mobile  erosion  flume  designed 
to  measure  cohesive  sediment  erosion  under  combined  wave- 
current  shear-stress  conditions.  Devices  have  been  previously 
developed  and  demonstrated  for  measuring  cohesive  sediment 
erosion  and  the  variation  with  depth  under  unidirectional,  steady 
current.  This  project  involves  developing  and  validating  a  similar 
flume  for  combined  current-wave,  unsteady  flow  conditions.  The 
flume  will  be  capable  of  measuring  erosion  variation  with  depth 
below  the  sediment/water  interface  so  that  the  well-documented 
increase  in  resistance  to  erosion  with  consolidation  is  measured. 

Achieving  this  technical  objective  will  address  the  need  for  cost-effective, 
accurate  methods  for  evaluating  the  potential  for  resuspension  of  buried 
contaminated  sediments  that  can  be  incorporated  into  risk  assessment 
framework. 


l 
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2  Background 

Aquatic  cohesive  sediments  are  often  the  final  receptor  of  contaminants 
released  by  Department  of  Defense  (DoD)  activities.  These  sediments  have 
accumulated  at  numerous  DoD  sites.  Understanding  the  stability  and 
erosion  of  these  fine  sediments  is  critical  to  assessing  contaminant  fate  and 
evaluating  remediation  options.  In  estuaries,  bays,  harbors,  and  channels, 
fine  sediments  can  be  mobilized  and  transported  by  the  combined  action  of 
both  relatively  steady  currents  (forced  by  river  flow,  tides,  and  wind)  and  to 
an  even  greater  extent  by  the  larger  unsteady  bottom  shear  stresses 
associated  with  higher-frequency  oscillatory  wave  motions  (either  wind-  or 
vessel-generated  waves).  Contaminated  site  remediation  is  very  expensive. 
Optimizing  impact  of  remediation  decisions  requires  defining  conditions 
under  which  buried  contaminated  sediment  may  erode. 

Cohesive  sediment  erosion 

The  majority  of  contaminated  sediments  are  a  mixture  of  sand,  silt,  and 
clay.  These  sediments  are  often  referred  to  as  cohesive  because,  unlike 
pure  sand,  bonding  between  particles  influences  sediment  behavior.  The 
complexity  of  cohesive  sediment  erosion  is  well  documented  (Gailani  et  al. 
2001;  Jepsen  et  al.  1997;  Jepsen  et  al.  2002;  McNeil  et  al.  1996;  Roberts  et 
al.  1998).  The  erosion  characteristics  of  these  sediments  cannot  be 
described  using  the  well-established  methods  developed  for  sandy 
sediments,  where  erosion  is  primarily  a  function  of  grain  size  distribution 
and  applied  shear  stress.  Cohesive  sediment  erosion  rates  will  be  affected 
not  only  by  grain  size  but  also  by  bulk  density,  mineralogy,  pore  water 
chemistry,  organic  content,  the  presence  of  gas  bubbles,  and  other  factors. 
It  is  qualitatively  understood  that,  depending  on  the  conditions,  small 
variations  in  one  or  more  of  these  bulk  properties  may  have  an  order  of 
magnitude  or  more  effect  on  erosion  rates.  However,  insufficient  data  are 
available  to  quantify  the  effects  of  these  properties  on  erosion  rates.  At 
present,  when  attempting  to  quantify  erosion  rates  and  conditions  for  a 
contaminated  cohesive  sediment  bed,  site-specific  sediment  erosion 
measurements  are  required  to  parameterize  erosion  equations.  In 
addition,  because  bulk  density  will  affect  erosion  rates,  these  rates  will 
change  with  depth  below  the  sediment-water  interface  even  for  well-mixed 
sediment,  adding  an  additional  level  of  complexity  to  erosion  testing.  Site- 
specific  measurement  methodology  exists  for  erosion  under  steady 
conditions;  however,  no  validated  devices  or  methods  exist  for  measuring 
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cohesive  sediment  erosion  (and  erosion  variation  with  depth)  under 
unsteady  conditions  similar  to  wave-current  environments. 

Measuring  cohesive  sediment  erosion 

The  U.S.  Army  Engineer  Research  and  Development  Center  (ERDC)  and 
Sandia  National  Laboratories  (SNL)  have  developed  two  portable 
field/laboratory  flumes  to  measure  cohesive  sediment  erosion  rate  variation 
with  depth  for  steady  flow  conditions.  The  first  flume  is  a  variation  of  the 
University  of  California-Santa  Barbara  (UCSB)  SEDflume  (McNeil  et  al. 
1996)  for  unidirectional,  steady,  high  shear-stress  flow.  Prior  to  SEDflume 
development,  only  surficial  sediment  erosion  could  be  reliably  measured 
using  portable  erosion  devices  (Tsai  and  Lick  1986;  Maa  et  al.  1993; 
Parchure  and  Mehta  1985).  Erosion  characteristics  of  underlying 
contaminated  sediments  measured  using  SEDflume  are  often  quite  different 
than  the  surficial  sediment  due  to  the  effects  of  consolidation,  biological 
activity,  and  other  processes.  A  second  flume,  the  ASSET  flume  (Jepsen  et 
al.  2005;  Roberts  et  al.  2003)  that  includes  all  the  capabilities  of  SEDflume 
as  well  as  the  capability  to  separate  bedload  and  suspended  load  transport 
of  the  eroded  material,  has  also  been  developed.  The  flumes  are  designed 
for  cost-effective  assessment  of  cohesive  sediment  erosion  rates.  Flume 
erosion  tests  replicate  the  applied  shear  stress  at  the  sediment  bed  (as 
opposed  to  flow  velocity)  and  relate  erosion  to  applied  stress  and  sediment 
bulk  properties.  One  major  benefit  of  replicating  shear  stress  instead  of 
velocity  regime  is  flume  size.  Shear  stress  can  be  replicated  in  a  relatively 
small  flume  and  with  a  small  sediment  sample  size.  Flumes  that  replicate 
flow  regime  must  be  significantly  larger,  require  a  larger  sediment  sample, 
and  are  generally  not  field  deployable. 

Cohesive  sediment  erosion  algorithms 

Algorithms  based  on  relationships  between  erosion  rate,  applied  shear 
stress,  and  sediment  bulk  properties  ultimately  allow  for  the  predictive 
modeling  of  bottom  sediment  fate  and  transport  and  sediment  bed  change. 
Extensive  SEDflume  application  was  used  to  develop  algorithms  for 
cohesive  sediment  erosion  variation  with  depth  under  unidirectional  flow 
(Roberts  et  al.  1998).  These  algorithms  have  been  demonstrated  to  be 
applicable  to  most  natural  cohesive  sediments.  These  algorithms  require 
site-specific  parameterization  developed  through  SEDflume  analysis  on  the 
sediments  of  interest.  Therefore,  flume  portability  to  the  field  is  considered 
essential  to  quantifying  erosion  rates  for  existing  sediment  beds  because 
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transporting  collected  sediment  cores  for  site-specific  testing  will  result  in 
an  excessive  disturbance  that  may  alter  erosion  behavior.  SEDflume  and 
ASSET  flume  permit  on-site  assessment  of  erosion  variation  with  depth  for 
minimally  disturbed  sediment  cores  extracted  from  the  contaminated  bed 
and  have  been  repeatedly  demonstrated  to  accurately  quantify  high-shear, 
cohesive  sediment  erosion  rate  variation  with  depth  through  numerous 
laboratory  and  field  applications  (Jepsen  et  al.  1997;  McNeil  et  al.  1996; 
Roberts  et  al.  1998).  Flume  data  are  used  to  develop  predictive  methods  for 
analyzing  the  stability  and/or  erosion  potential  of  cohesive  sediment  beds, 
including  stratified  beds.  These  predictions  are  used  in  risk  assessment,  site 
management,  and  assessment  of  remediation  options  for  contaminated 
sediment  sites  (Borrowman  et  al.  2006).  The  Environmental  Protection 
Agency  (EPA),  for  example,  utilizes  SEDflume  at  most  Superfund  sites 
where  buried  contaminated  sediment  deposits  are  an  issue.  The  flume  is 
also  commonly  used  in  developing  dredged  material  management  plans  and 
assessing  the  impact  of  cohesive  sediment  deposits  on  surrounding 
ecosystems. 

Cohesive  sediment  erosion  under  unsteady  flow  conditions 

Unidirectional,  flow-induced,  shear-stress  conditions  in  the  SEDflume  and 
ASSET  flume  are  a  reasonable  representation  of  shear-stress  conditions  in 
rivers  or  in  deep  water  where  wave  influence  on  bottom  shear  stress  is  not 
significant.  However,  many  contaminated  sediment  sites  are  in  the  shallow 
waters  of  lakes,  reservoirs,  harbors,  large  rivers,  estuaries,  or  coastal  regions 
where  wave  action  influences  near-bottom  shear  stresses.  SEDflume  and 
ASSET  flume  erosion  rates  measured  under  steady,  unidirectional  flow  are 
not  applicable  to  these  environments.  Therefore,  additional  flumes  are 
required  that  represent  the  time-varying  conditions  experienced  in  wave- 
influenced  environments.  The  ERDC  and  SNL  recently  developed  a 
prototype  flume  for  measuring  erosion  under  unidirectional  flow  with  a 
superimposed  oscillatory  motion;  these  conditions  represent  shear 
conditions  under  a  combined  wave-current  environment.  The  new  flume  is 
the  previously  defined  SEAWOLF  (Jepsen  et  al.  2002;  Kearney  et  al.  2008). 
SEAWOLF  is  based  on  the  proven  SEDflume  and  ASSET  flume  designs. 

Like  its  predecessors,  SEAWOLF  mimics  the  shear-stress  time  history 
induced  by  given  site  conditions,  as  opposed  to  mimicking  the  actual 
hydrodynamics  of  the  field,  and  will  be  a  cost-effective,  field-deployable 
method  for  measuring  erosion  rate  variation  with  depth  for  minimally 
disturbed  core  samples.  Algorithms  based  on  sediment  parameterization 
achieved  by  SEAWOLF  flume  testing  will  be  developed  in  this  research 
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project  so  that  contaminated  sediment  fate  and  transport  can  be 
predicatively  modeled  in  combined  wave-current  environments.  To  this 
end,  the  new  flume  requires  validation  of  applicability,  similar  to  the 
multiyear  effort  required  to  validate  SEDflume. 

Description  of  SEAWOLF  flume 

SEAWOLF  features  a  2  m  long,  straight,  clear  polycarbonate  flume 
channel  (Figure  1).  The  channel  has  a  false  bottom  at  the  center  where  a  10 
cm  diameter  sediment  core  extracted  directly  from  the  field  site  (or 
created  in  the  laboratory)  is  placed  (Figure  2).  Sediment  cores  can  be  up  to 
1  m  in  length  below  the  initial  sediment  bed  surface.  As  with  SEDflume, 
the  sediment  core  is  eroded  within  the  flume  channel  test  section  to 
quantify  erosion  rate.  Sediment  traps  at  each  end  of  the  flume  channel 
remove  sediments  from  the  system  such  that  core  erosion  in  the  test 
section  is  not  influenced  by  sediment-laden  water  or  bedload  from 
previously  eroded  material. 


Figure  1.  SEAWOLF  flume  schematic.  Three-tank  system  and  flume  that  drives  the  unidirectional  flow  via  a  head 
difference  (Ah)  between  Tank  A  and  Tank  C.  Tank  B  is  overflow  from  Tanks  A  and  C,  which  is  sent  back  to  Tank 

A  via  the  pump. 
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Figure  2.  Bottom  view  of  circular  insert  hole  in  base  plate  where 
sample  core  is  inserted.  Rubber  seals  have  been  waterproofed  with 
silicon  and  secured  with  nuts  and  hex  bolts. 


Tanks  A  and  C  are  connected  to  opposite  ends  of  the  flume  channel  test 
section.  Water  is  pumped  from  Tank  B  to  Tank  A  to  maintain  the  desired 
head  in  Tank  A.  The  head  in  Tank  A  is  greater  than  the  head  in  Tank  C. 
This  head  difference,  Ah,  which  drives  the  unidirectional  flow  in  the 
SEAWOLF  test  section,  can  be  adjusted  between  each  erosion  test  to 
represent  user-specified  mean  current  velocities.  Both  Tank  A  and  Tank  C 
overflow  into  Tank  B  to  maintain  constant  Ah  during  an  erosion  test.  A 
magnetic  flow  meter  attached  directly  to  one  end  of  the  channel  provides 
real-time  measurements  of  flow  conditions. 

Two  pistons  (Figure  3)  working  in  tandem  (Figure  1)  generate  an 
oscillatory  flow  over  the  test  section  that  represents  wave-induced  orbital 
velocities  at  the  sediment  bed.  A  computer,  stepper  motor,  and  linear  ball- 
screw  arrangement  control  the  piston  strokes  that  govern  the  velocity  time 
history  and  period  of  the  oscillatory  flow.  In  addition,  valves  at  each  end  of 
the  channel  connecting  to  Tank  A  or  Tank  C  are  used  to  control  both  the 
unidirectional  flow  rates  and  the  backflow  into  the  tanks  from  the 
oscillatory  flow.  Piston  movement  and  head  difference  are  pre¬ 
programmed  by  the  operator  before  each  erosion  test. 
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Figure  3.  Dual  pistons,  which  create  the  oscillatory  flow  in  the  flume. 
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The  rectangular,  cross-section  test  channel  is  2  m  long  and  20  mm  high  x 
105.4  mm  wide  (aspect  ratio  =  5.5:1).  All  PIV  measurements  are  made  in  a 
streamwise/ wall-normal  plane,  which  is  located  at  the  spanwise  centerline 
at  a  streamwise  location  of  1  m,  or  52.5  channel  heights,  from  the  channel 
ends.  BK-7  optical  flats  of  50.8  mm  diameter  are  mounted  flush  with  the 
test-section  floor  and  ceiling  at  the  measurement  location  to  allow  for 
passage  of  the  Nd:YAG  laser  sheet  used  in  the  PIV  measurements.  The 
test-section  ceiling,  floor,  and  sidewalls  are  fabricated  from  12.7  mm  thick 
acrylic.  A  BK-7  optical  flat  is  mounted  in  one  of  the  sidewalls  for  PIV 
imaging  at  right  angles  to  the  laser  sheet. 

The  test  channel  was  sufficiently  long  to  provide  fully  developed  mean 
velocity  profiles  under  steady  flow  conditions,  as  confirmed  by  PIV- 
measured  velocity  profiles  that  conformed  to  the  well-established  law  of 
the  wall,  and  wall  shear  stresses  that  were  generally  within  10%  of  the 
Blasius  correlation  (White  1991)  for  fully  developed  turbulent  wall- 
bounded  flow.  The  2  m  test-section  length  was  limited  by  space 
constraints  and  was  short  enough  that  end  effects  were  present,  as  fluid 
that  began  the  forcing  cycle  at  the  measurement  station  left  and  re-entered 
the  test  section  during  the  course  of  a  forcing  cycle. 
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3  Materials  and  Methods 

Based  upon  the  very  interesting  results  from  the  hydrodynamic  data 
collection  activities,  and  with  the  SERDP  review  panel’s  blessing,  the  path 
forward  changed  in  2008  to  focus  exclusively  on  hydrodynamic  calibration 
and  validation.  Therefore,  the  objective  evolved  to  include  a  series  of 
experiments  to  validate  the  applicability  of  the  prototype  SEAWOLF  flume 
to  replicate  shear-stress  time  histories  in  naturally  occurring  wave-current 
conditions.  The  change  in  path  forward  has  necessitated  a  redevelopment 
of  the  milestones  for  this  research  project.  The  new  tasks  and  subtasks  are 
outlined  below. 

•  Task  1  -  Determine  shear-stress  time  histories  and  turbulent  structure 
within  the  enclosed  SEAWOLF  channel  with  smooth  and  rough  walls 
under  modulated  flow 

1.1  Assess  the  need  to  modify  SEAWOLF  Flume  Channel  to  contain 
all  oscillating  flow  and  modify  as  necessary 

1.2  PIV  measurements  in  smooth-walled  channel  flows 

1.3  PIV  measurements  in  various  rough-walled  channel  flows 

•  Task  2  -  Create  data  base  of  shear-stress  time  histories  by 
demonstrating  and  applying  newly  developed  optical  sensing 
technology  within  SEAWOLF 

2.1  Identify  and  test  capabilities  of  an  optical  shear-stress  sensor 
within  SEAWOLF 

2.2  Implementation  of  optical  shear-stress  sensor  in  smooth-walled 
channel  flow 

2.3  Optical  shear-stress  measurements  in  rough-walled  channel  flow 

2.4  Generate  comprehensive  database  of  shear-stress  time  histories  in 
SEAWOLF  using  optical  shear-sensor  technology  measurements 

•  Task  3  -  Demonstrate  wave-current  bottom  shear-stress  time  histories 
representative  of  nature  in  SEAWOLF  by  manipulating  flume  and 
modulated  flow  conditions 

3.1  Assessment  of  existing  knowledge  of  shear-stress  time  histories  in 
naturally  occurring  bottom  boundary  layers 
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3.2  Configure  SEAWOLF  Flume  for  replication  of  wave-current 
boundary  layer  shear  stress 

3.3  Replicate  naturally  occurring  wave-current  conditions  in 
SEAWOLF  Flume 

•  Task  4  -  Investigate  uncertainty  in  shear-stress  measurement 

techniques 

4.1  Assess  uncertainty  in  PIV  shear-stress  measurements 

4.2  Assess  uncertainty  of  optical  sensor  shear-stress  measurements 

These  tasks  outline  the  specific  objective  of  this  research,  which  is  to  focus 
on  quantifying  and  controlling  time-varying  bottom  shear  stress  within  the 
enclosed  channel  of  the  SEAWOLF  for  combined  unidirectional  and 
oscillatory  flow  (modulated  flow).  These  combined  flows  are  representative 
of  near-bottom  hydrodynamic  conditions  in  a  wave-current  environment 
and  are  therefore  representative  of  shear  stress  applied  to  the  sediment  bed. 
Velocity  profile  and  turbulence  length  scale  have  not  been  satisfactorily 
measured  in  an  enclosed  channel  for  the  purpose  of  quantifying  wall  (bed) 
shear  stress  in  modulated  flow.  Therefore,  it  is  imperative  to  measure 
phase-resolved  turbulent  flow  structure  and  wall  shear  stress  within  a 
confined  flume  with  modulated  flow.  Once  measurement  techniques  have 
been  demonstrated,  the  next  step  is  to  control  the  turbulent  flow  structure 
to  construct  a  desired  shear-stress  time  history  that  is  relevant  to  nature  by 
adjusting  bottom  roughness,  flume  structure,  and  flow  parameters. 

Task  1-Determine  shear-stress  time  histories  and  turbulent  structure 
within  the  enclosed  SEAWOLF  channel  with  smooth  and  rough  walls 
under  modulated  flow 

Erosion  is  generally  quantified  as  a  function  of  bottom  shear  stress.  Shear 
stress  is  a  function  of  the  near-bottom  velocity  profile  (not  mean  velocity) 
and  sediment  bed  roughness.  The  complex,  time-varying  bottom  shear 
stress  in  SEAWOLF  had  not  been  extensively  measured  or  tested  prior  to 
Strategic  environmental  Research  and  Development  Program  (SERDP) 
funding.  To  accurately  quantify  wave-current  induced  erosion,  the 
SEAWOLF  bottom  shear-stress  time  series  within  the  flume  must  reflect 
the  time  series  of  shear  stress  due  to  a  specified  combined  current-wave 
condition  that  occurs  in  nature.  Currently,  theory  cannot  describe  time 
series  shear  stress  within  an  enclosed  channel  subject  to  oscillating  flow 
(theory  can  be  and  is  used  to  relate  flowrate  to  shear  stress  in  SEDflume 
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and  the  ASSET  flume— no  measurements  were  deemed  necessary  in  these 
steady  flow  flumes).  Shear-stress  time  history  in  the  flume  channel  must 
be  quantified  through  direct  measurement  of  near-bottom  velocity  profile, 
leading  to  a  calculation  of  bottom  shear  stresses.  These  measurements 
permit  control  of  piston  movement  and  head  difference  to  replicate 
bottom  shear-stress  time  histories  generated  during  specified  wave- 
current  conditions. 

Particle  image  velocimetry  (PIV)  operating  principles 

The  PIV  technique  is  described  in  picture  form  in  Figure  4.  With  PIV,  the 
flow  of  interest  is  seeded  with  particles  that  are  both  small  and  light 
enough  to  faithfully  track  the  flow.  The  particles  are  illuminated  using  a 
planar  sheet  of  pulsed  laser  light.  The  laser  pulses  must  be  short  enough  to 
effectively  freeze  the  particle  field,  such  that  particle  motion  during  the 
laser  illumination  is  negligible.  For  this  purpose,  pulsed  Nd:YAG  lasers, 
that  provide  visible  532  nm  (green)  light  with  a  short,  10  ns  pulse 
duration,  are  typically  used.  Simple  cylindrical  and  spherical  lens  elements 
are  typically  used  to  re-form  the  laser-light  output  into  a  thin  sheet  of 
light,  which  illuminates  the  tracer  particles  in  the  flow.  For  PIV 
measurement,  two  laser  pulses  are  delivered  in  succession  with  a  known, 
prescribed  time  delay,  At,  between  the  laser  pulses.  An  interline  transfer 
CCD  camera  is  then  used  to  record  two  digital  particle  images  in  rapid 
succession,  with  exposures  which  are  timed  such  that  each  imag e  frames 
one  of  the  laser  pulses.  The  particle-image  pairs  contain  information  about 
an  instantaneous  snapshot  of  the  velocity  field  in  the  plane  of  the  laser 
sheet.  The  digital  images  are  broken  down  into  16-32  pixel  square 
interrogation  windows  and  the  average  displacement,  Ax,  of  the  particles 
in  each  window  during  the  time  At  between  the  laser  pulses  is  determined 
using  statistical  cross-correlation  techniques.  The  velocity  vector  in  each 
interrogation  spot  is  then  determined  from 


11  =  Ax/ A/ 


(1) 
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Figure  4.  Conceptual  depiction  of  the  PIV  system. 


Laser  Illumination 


/ 

cross 

correlation 

•  • 

•  • 

•  • 

•  • 

• 

• 

v=Ax/At  t=t0  t=t0+At 

Analysis  of  Local  Interrogation  Spots 


Since  both  the  horizontal  and  vertical  particle  displacements  of  the  tracer 
particles  can  be  determined  from  the  two-dimensional  (2D)  particle 
images,  u  is  then  a  vector  quantity  containing  the  two  components  of  the 
velocity  field,  u  and  v  (horizontal  and  vertical),  in  the  plane  of  the  laser 
sheet.  When  this  process  is  repeated  at  each  interrogation  spot  in  the 
images,  2D  vector  maps,  such  as  those  shown  in  Figure  4  are  obtained. 
These  maps  reveal  the  instantaneous  structure  of  a  complex  turbulent 
flow.  Many  such  maps  can  then  be  statistically  analyzed  to  extract  the 
mean  velocity  field  and  the  statistical  fluctuations  of  the  velocity  which 
define  the  turbulent  Reynolds  stress  tensor. 

PIV  optical  system  and  tracer  particles 

A  schematic  of  the  PIV  optical  system  is  presented  alongside  the  drawing  of 
the  SEAWOLF  channel  in  Figure  5.  A  dual-cavity,  frequency-doubled 
Nd:YAG  laser  (Continuum  MiniLite)  provides  Q-switched  laser  pulses  of  8 
ns  duration,  with  energies  of  -25  mJ/pulse  at  a  wavelength  of  532  nm. 
Particle  images  were  captured  on  an  interline-transfer  CCD  camera  (Kodak 
ES  1.0)  with  1024  x  1024  pixel  format  and  subsequently  stored  for  cross¬ 
correlation  postprocessing.  PIV  experiments  at  two  magnification  settings 
were  conducted:  (1)  a  low-magnification  (0.33:1)  setting,  attained  using  a 


ERDC/CHL  TR-17-15 


12 


105  mm  Nikkor  lens  at//8  aperture,  was  used  to  obtain  PIV  recordings  with 
a  field  of  view  that  spanned  the  full  height  of  the  channel  and  (2)  a  high- 
magnification  (-3.45:1)  setting,  using  the  same  105  mm  lens  with  a  2x 
teleconverter  and  extension  rings,  to  obtain  a  field  of  view  which  was 
nominally  2.5  mm  square.  The  low-magnification  PIV  experiments  allowed 
viewing  of  the  large-scale  behavior  of  the  unsteady  flow  field  while  the  high- 
magnification  data  sets  permit  resolution  of  the  viscous  sublayer  (described 
in  detail  in  subsequent  sections)  and  obtain  a  direct  measurement  of  the 
wall  shear-stress  modulation  throughout  the  cycle. 


Figure  5.  Schematic  of  the  SEAWOLF  test  facility 
and  PIV  Instrumentation. 


The  flow  was  seeded  with  1.03  g/cc,  5-//m-diameter  polyamide  particles, 
which  were  nearly  neutrally  buoyant  in  water.  The  SEAWOLF  facility  was 
further  contaminated  by  residual  sediment  particles  from  earlier  testing, 
which  were  estimated  to  be  as  large  as  35  /an,  and  could  not  be  completely 
removed.  The  flow-tracking  ability  of  both  seed  and  contaminant  particles 
was  estimated  from  the  particle  Stokes  number: 

,2 

pn  a„ 

Sk  =  ^-^fc  (2) 

18  n 


where  pp  is  the  particle  density,  dp  is  the  particle  diameter,  //  is  the 
dynamic  viscosity,  and/c  is  a  characteristic  frequency  of  the  flowfield  that 
can  be  estimated  from  the  rates  of  strain  present  in  the  flow.  For  a  worst- 
case  estimate,  fc  was  computed  using  the  viscous  timescale,  t+  ,  where 

Ur  is  the  friction  velocity,  and  vthe  kinematic  viscosity,  as  t+  provides  an 
estimate  of  the  strain  rate  in  the  high-velocity-gradient  region  at  the  wall. 
The  polyamide  particles,  as  expected,  yielded  Sk  values  of  order  io~3, 
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compared  to  values  of  o.i  to  0.7  for  the  sediment  particles.  These  worst- 
case  values  were  small  enough  to  guarantee  faithful  particle  tracking  in 
most  cases,  with  the  exception  of  sediment  contaminants  in  the  very  near¬ 
wall  region  and  only  at  the  high  end  of  the  shear  stresses  present  in  this 
study.  As  discussed  in  section  IV,  particle  tracking  does  not  impact  the 
estimates  of  wall  shear  stress,  which  were  validated  using  a  steady 
turbulent  channel  flow  as  a  test  case. 

Representative  high-  and  low-magnification  particle  images  are  shown  in 
Figure  6.  The  low-magnification  images  are  of  excellent  quality  and 
permitted  extraction  of  high-resolution  PIV  results.  Good  particle-image 
density  was  observed  in  the  high-magnification  realizations,  but  the 
quality  of  the  high-magnification  results  was  degraded  somewhat  as  the 
diffraction-limited  spot  size  becomes  comparable  or  better  than  the 
particle  diameters,  resulting  in  a  greater  disparity  in  particle-image  size 
between  seed  and  contaminant  particles.  No  processing  of  the  particle 
images  was  performed  to  differentiate  between  small  and  large  particle 
images,  and  velocity  data  were  necessarily  obtained  from  both  seed  and 
contaminant  particles. 

Figure  6.  Representative  particle  images  at  low-  and  high-magnification  settings. 
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Laser  timing  for  unsteady  flow  measurements 

The  large  variation  in  bulk  velocity  through  the  cycle  required  a  continuous 
varying  of  the  time  delay  between  the  Nd:YAG  laser  pulses  used  for  PIV.  A 
timing  diagram  for  the  PIV  measurements  is  provided  in  Figure  7  to 
illustrate  this  procedure.  All  timing  pulses  were  generated  using  counter/ 
timer  chips  provided  with  commercial  data-acquisition  cards,  controlled 
using  a  custom  LABVIEW  program.  A  fixed  number  of  40  preprogrammed 
laser-pulse  delays  were  manually  entered  into  a  list  in  the  LABVIEW 
software,  which  was  coded  to  provide  a  set  of  40  triggers  per  flow  cycle,  and 
then  there  was  a  wait  until  the  next  cycle  was  initiated  by  the  SEAWOLF 
pistons.  The  start  of  the  flow-forcing  cycle  was  indicated  by  the  transistor- 
transistor  logic  (TTL)  pulse  (uppermost  in  Figure  5)  from  an  optical  switch, 
which  was  actuated  at  the  end  of  travel  of  one  of  the  pistons  when  the  path 
between  its  light  source  and  photodiode  was  broken  by  a  metal  tab  that 
traveled  with  the  piston  shaft.  Trigger  pulses  for  the  CCD  camera  and  laser 
A  (second  and  third  pulse  trains  in  Figure  7)  were  provided  in  uniform  r/40 
increments,  where  T  is  the  period  of  the  flow  forcing  cycle.  The  delay 
between  laser  A  and  laser  B  was  then  varied  (fourth  pulse  train  in  Figure  7) 
using  the  preprogrammed  list  to  the  LABVIEW  code,  which  cycled 
automatically  through  the  sequence  and  then  stood  by  until  the  beginning  of 
the  next  forcing  cycle  was  triggered  by  the  optical  switch.  The  trigger  to  the 
CCD  camera  was  also  used  to  acquire  data  from  a  volume  flow  meter  and  a 
thermocouple  located  near  the  entrance  to  the  test  section  simultaneously 
with  the  PIV  recordings.  The  variable  delay  between  laser  pulses  was  set  to 
obtain  a  maximum  streamwise  displacement  of  nominally  10-15  pixels  in 
the  particle-image  pairs,  which  yields  a  dynamic  velocity  measurement 
range  of  100-150:1  for  the  streamwise  component  and  10-20:1  for  the 
vertical  component,  based  on  a  commonly  assumed  (Raffel  et  al.  2002)  ±0.1 
pixel  uncertainty  in  the  location  of  the  cross-correlation  peak. 


Figure  7.  Timing  diagram  for  PIV  measurements  in  oscillating  flow. 
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PIV  image  processing  and  spatial  resolution 

Our  particle  images  were  interrogated  using  cross-correlation  algorithms 
implemented  in  the  DaVis  commercial  software  package  (DaVis  2006).  The 
procedure  was  initiated  using  a  coarse  64  x  64  pixel  interrogation-window 
size  with  no  window  offset,  which  was  required  since  the  direction  of  the 
streamwise  velocity  was  not  known  a  priori.  The  coarse  field  was  liberally 
interpolated  and  smoothed  to  provide  an  estimate  for  window  displacement 
in  subsequent  passes  with  finer  resolution.  The  final  PIV  result  was 
obtained  using  16  x  16  pixel  windows  with  50%  spatial  overlap.  Outliers 
were  removed  from  the  vector  fields,  and  holes  in  the  final  vector  estimates 
were  conservatively  interpolated.  The  spatial  resolution  in  the  velocity 
vector  data  was  160  /an  and  19.6  /an  for  the  low-  and  high-magnification 
measurements,  respectively.  When  scaled  to  inner  variables,  the 
corresponding  in-plane  spatial  resolutions  range  from  1.7  to  13.1  viscous 
wall  units  for  the  low-magnification  data  and  0.2  to  1.6  wall  units  at  high- 
magnification,  for  a  measured  range  of  wall  shear  stresses  of  0.1  to  6  Pa.  The 
out-of-plane  resolution  in  the  data  was  estimated  at  2  to  20  wall  units, 
based  on  an  assumed  200  /an  overlap  of  the  laser  sheets. 


Supporting  theory  for  PIV-derived  wall  shear-stress  measurements 

A  key  objective  of  this  effort  is  the  determination  of  wall  shear  stress  by 
direct  measurement  of  the  gradient  of  the  phase-resolved  velocity  profile 
at  the  wall,  using  the  well-known  constitutive  relation: 


T 


W 


(3) 


where  rw  is  the  wall  shear  stress,  //  is  the  dynamic  viscosity,  and  du/dy  is  the 
slope  of  the  mean  velocity  profile  at  the  wall,  which  is  to  be  estimated  from 
high-magnification  PIV  data.  For  this  procedure  to  be  valid,  a  sound 
justification  is  sought  for  a  linear  fit  to  the  PIV-measured  velocity  data  in 
some  region  close  to  the  wall.  The  following  analysis  is  a  review  of  well- 
known  dimensional  scaling  laws  for  flow  in  the  innermost  portions  of  a 
boundary  layer,  very  close  to  the  wall,  as  presented  by  Hinze  (1987) 1.  This 
analysis  justifies  the  use  of  a  linear  representation  of  the  velocity  profile  in 
the  near-wall  region,  which  makes  PIV  shear-stress  determination  possible. 


1  See  section  7-5  of  Hinze  (1987)  for  a  discussion  of  scaling  in  the  inner  portion  of  the  boundary  layer. 
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Very  close  to  the  wall,  the  dynamics  of  the  boundary  layer  flow  are 
dominated  by  the  momentum  transfer  at  the  wall,  through  the  shear 
stress,  given  by  Equation  3.  If  Equation  3  is  divided  by  the  density,  p,  a 
velocity  scale  that  is  dependent  on  the  shear  stress  is  obtained,  which 
drives  the  flow  in  the  near  wall  region: 


tw  _ydu 
P  dy 


(4) 


where  vis  the  kinematic  viscosity  and  uT  has  units  of  velocity  and  is 
commonly  referred  to  as  the  friction  velocity,  which  is  the  appropriate 
velocity  scale  for  the  near-wall  flow.  Dimensional  analysis  of  Equations  2 
and  3  reveals  that  the  near- wall  velocity  profile  u  =  u(uT,v,y),  which  can 
be  reduced  to  a  functional  relationship  between  dimensionless  variables 
using  the  n  theorem: 


u 

Ur 
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uTy 


{  V 


(5) 


or 


u+  =  fly) 


(6) 


Equation  (6)  is  a  functional  form  or  scaling  law,  which  holds  in  the  so- 
called  inner  region  of  a  boundary  layer,  very  far  from  the  boundary  layer 
edge,  channel  centerline,  or  other  boundary,  so  that  the  shear  stress  at  the 
wall  solely  governs  the  flow  dynamics.  As  a  result,  u+  and  y+  are  referred  to 
as  inner  variables. 

The  form  of  the  function/ (y+)  will  dictate  the  equation  that  is  needed  to  fit 
the  PIV  data  in  the  near-wall  region.  In  the  immediate  vicinity  of  the  wall, 
the  no-slip  condition  requires  that  the  velocity  be  near  zero.  Transport  of 
momentum  is  dominated  by  the  molecular  collisions  represented  by  the 
viscosity  so  that  the  velocity  profile  is  entirely  determined  by  Equation  3, 
which  is  integrated  from  the  wall  to  a  near-wall  position,  y,  to  obtain  the 
following  expression  for  the  near-wall  velocity  profile: 
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which  can  be  rearranged  by  dividing  by  urto  obtain 


u+  =  y+  for  y+  <  5 


(8) 


Equations  7  and  8  are  valid  in  the  very  low-momentum  region 
immediately  adjacent  to  the  wall,  where  viscous  diffusion  is  the  dominant 
force  affecting  the  velocity  profile.  Experimentally,  Equation  6b  has  been 
shown  to  be  valid  for  only  a  very  thin  region  between  the  wall  and  y+  =  5, 
beyond  which  turbulent  motion  begins  to  impact  the  transport  in  addition 
to  purely  viscous  forces.  This  relationship  reveals  that  the  velocity  profile 
in  this  viscous  sublayer  is  linear  and  that  a  linear  fit  to  PIV  data  measured 
very  near  the  wall  is  justified  for  the  determination  of  du/dy  and  the  wall 
shear  stress  through  Equation  3. 

Proceeding  farther  from  the  wall,  turbulent  mixing  just  begins  to  impact 
the  transport  of  momentum  in  the  near-wall  region  and  becomes 
comparable  or  even  greater  than  the  molecular  viscous  forces.  The  impact 
of  this  turbulent  mixing  manifests  itself  as  an  apparent  time-averaged 
shear  stress,  which  can  be  expressed  in  a  form  similar  to  Equation  3  using 
a  turbulent  or  eddy  viscosity,  vt.  The  combined  impact  of  the  viscous  and 
turbulent  shear  forces  can  then  be  written  as 


r 

P 


+  T 


turbulent 


P 


(v+vt) 


du 

dy 


(9) 


Experiments  show  that  by  a  distance  of  approximately  y+  =  30,  the 
turbulent  mixing  dominates  and  vt»  v.  Based  on  physical  arguments 
(Akhavan  et  al.  1991)  vt  can  be  expressed  as 

vt=KUTy  (10) 


where  /cis  a  proportionality  constant.  If  Equation  10  is  substituted  into 
Equation  7  and  the  turbulent-mixing-dominated  case  is  considered,  with 
vt»  v,  for  the  region  for  y+  >  30,  the  following  is  obtained: 


r 

P 


2  du 

ur=K^y~r 

dy 


(id 


or  in  dimensionless  quantities, 
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*y 


du+ 

dy+ 


(12) 


Equation  12  can  be  integrated  to  obtain  the  following  functional  form  for 
the  mixing-dominated  region  of  the  inner  layer: 

u+  =  —  \ny+  +B  fory*  >  30  (13) 

K 

For  steady,  wall-bounded  turbulent  flows,  the  constants  have  been 
determined  experimentally  to  be  k  =  0.41  and  B  ~  5-6. 

Equations  7  and  13  represent  analytical  forms  for  turbulent  boundary 
layers.  Since  the  dimensionless  variables,  u+  and  y+,  in  these  expressions 
depend  on  the  value  of  the  wall  shear  stress,  these  equations  can  be  used  to 
fit  the  near-wall  PIV  data  for  the  determination  of  the  wall  shear  stress.  The 
linear  expression  for  the  viscous  sublayer  in  Equation  8  is  valid  for  steady 
and  unsteady,  laminar  and  turbulent  flows.  Furthermore,  Equation  6b  has 
no  empirically  determined  constants  and  is  thus  a  general  way  of 
determining  shear  stress  from  measured  velocity  profiles.  The  drawback  to 
Equation  8is  that  the  viscous  sublayer  is  often  very  thin,  and  velocity 
measurements  very  close  to  the  wall  are  extremely  challenging.  The 
logarithmic  region  for  y+  >  30  is  simpler  to  access  experimentally,  and  fits 
to  the  form  of  Equation  10  are  often  used  for  the  measurement  of  wall  shear 
stress  in  steady  turbulent  boundary  layers  (Clauser  1956;  Kendall  and 
Koochesfahani  2006).  The  derivation  of  Equation  10  presented  here  is  for 
steady  turbulent  flows  only,  but  Akhavan  et  al.  (1991)  has  shown  that  the 
logarithmic  form  is  also  valid  for  unsteady  turbulent  boundary  layers  as 
well.  The  two  main  drawbacks  to  logarithmic  fitting  in  unsteady  flows  are 
(1)  the  log  law  only  applies  under  turbulent  conditions,  while  the  oscillating 
and  pulsating  flows  considered  here  undergo  periods  of  laminar  and 
transitional  flow,  and  (2)  the  constants  /rand  B  in  the  log  law  are  not  well 
known  for  unsteady  flows,  and  they  can  change  throughout  the  forcing 
cycle,  further  making  Equation  10  difficult  to  apply  in  unsteady  flows. 

Fitting  of  the  viscous  sublayer  and  use  of  Equation  3  was  determined  to  be 
the  best  option  for  PIV-based  shear-stress  measurements  in  SEAWOLF,  as 
this  method  was  most  general  and  likely  to  be  the  most  accurate.  With  this 
in  mind,  the  high-magnification  PIV  instrument  previously  described  was 
constructed.  For  the  0-6  Pa  range  of  shear  stresses  observed  in  the 
SEAWOLF  water  flow,  the  thickness  of  the  viscous  sublayer  varies  from 
approximately  70  to  220  jum  so  that  any  PIV  measurement  for  wall-shear- 
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stress  determination  must  have  very  fine  spatial  resolution.  The  19.6  /an 
spatial  resolution  of  the  high-magnification  system  permits  3  to  11  points 
within  the  viscous  sublayer,  to  which  a  linear  fit  can  be  performed  for  the 
estimation  of  the  wall  velocity  gradient,  and  direct  measurement  of  the 
wall  shear  stress  through  Equation  3.  The  fidelity  of  this  approach  was 
then  verified  with  measurements  in  steady  turbulent  channel  flows,  as 
explained  in  the  Task  1.2  subsection  of  the  Results  and  Accomplishments 
section  of  this  report. 

Steady  flow  validation  of  PIV  wall  shear-stress  measurements 

A  key  objective  of  this  effort  is  the  determination  of  wall  shear  stress  by 
direct  measurement  of  the  gradient  of  the  phase-averaged  velocity  profile  at 
the  wall  using  the  high-magnification  PIV  approach  suggested  by  Oweis  and 
Ceccio  and  employed  in  their  work  (Oweis  et  al.  2010).  The  0.2  to  1.6 
viscous-wall-unit  resolution  of  the  high-magnification  PIV  data  enables 
probing  of  the  viscous  sublayer  for  the  Reynolds  numbers  of  interest  here. 
Wall  shear  stress  is  calculated  from  Equation  3  by  estimating  the  velocity 
gradient  from  the  slope  of  a  least-squares  line  fit  to  all  mean  velocity  data 
with  i/+  <  5.  The  effectiveness  of  this  procedure  was  investigated  using  data 
from  steady  turbulent  flows  established  in  the  SEAWOLF  test  section.  The 
shear-stress  data  were  then  compared  with  redundant  estimates  of  rw 
obtained  from  Clauser-type  fits  layers  (Clauser  1956;  Kendall  and 
Koochesfahani  2006)  to  the  velocity  data  in  the  logarithmic  layer  and  from 
the  empirical  shear-stress  correlation  due  to  Blasius  (White  1991): 

^-=0.0791^1/4  (14) 

pU 

where  p  is  the  density  and  U  is  the  mass-mean  or  bulk  velocity  in  the 
channel. 

These  steady-flow,  high-magnification  PIV  results  are  shown  in  Figure  8. 
Mean  velocity  profiles  were  computed  from  ensemble  and  streamwise 
averaging  of  500  statistically  independent  PIV  realizations,  as  shown  in 
the  uppermost  plot  in  Figure  8.  The  mean  velocity  data  were  normalized 
using  inner  length  scales  computed  from  fitting  of  the  log  layer.  At  all 
Reynolds  numbers,  the  velocity  profiles  display  good  collapse  to  the 
universal  curve,  with  some  scatter  observed  in  the  results  closer  to  the 
wall,  which  is  indicative  of  uncertainty  in  wall  location  and  small 
differences  in  the  viscous  scales  computed  using  shear-stress  values  from 
log-layer  and  viscous-layer  fitting. 
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Figure  8.  High-magnification  PIV  results  from  a 
steady  turbulent  boundary  layer:  (a) 
streamwise  velocity  profiles  expressed  in 
terms  of  inner  variables,  (b)  wall  shear-stress 
data  and  comparison  to  Blasius  correlation. 


Wall  shear-stress  data  from  viscous  and  log-layer  fits  are  provided  as  a 
function  of  Reynolds  number  in  the  lowermost  plot  in  Figure  8.  The 
experimental  results  are  plotted  alongside  empirical  values  computed 
from  Equation  14.  Shear-stress  estimates  from  linear  fits  to  the  velocity 
data  in  the  viscous  layer  are  within  -4%  to  -13%  of  both  Equation  6  and 
results  from  log-layer  fitting,  which  was  seen  as  reasonable  level  of 
accuracy  for  the  measurements  in  the  oscillating  channel  flow.  Note  that 
the  outstanding  agreement  between  the  log-layer-fit  results  and  Equation 
11  is  surely  fortuitous  and  that  none  of  the  three  shear-stress  estimates  in 
Figure  8  can  be  deemed  as  more  reliable  than  any  other.  It  is  the  ~io% 
agreement  among  all  three  shear-stress  estimates  that  gives  confidence  in 
the  approach  of  fitting  to  the  viscous-layer  data. 

Phase  averaging  procedures  for  velocity  and  shear-stress  results 

All  velocity  and  turbulence  quantities  are  presented  in  the  context  of  a 
phase-averaging  procedure  that  is  common  to  the  analysis  of  periodically 
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unsteady  turbulent  flows  (Akhavan  et  al.  1991;  Jensen  et  al.  1989;  Hino  et 
al.  1983).  In  this  analysis,  any  quantity,  q  (representing  a  velocity 
component,  shear  stress,  or  other  variable),  can  be  expressed  on  an 
instantaneous  basis  in  terms  of  a  phase-locked  average  (9)  and  an 
instantaneous  turbulent  fluctuation,  q' : 

q(x,T)  =  (q)(x,T)  +  q'(x,t)  (15) 


where  the  functional  dependences  on  position  vector,  x;  the  period- 
normalized  time,  t  =  t/T  (ranges  from  o  to  1  over  the  course  of  the  cycle); 
and  an  instantaneous  time,  t,  are  explicitly  indicated.  From  discretely 
sampled  data,  the  phase  average  simply  provides  a  representation  of  the 
average  cycle,  and  can  be  computed  from 


<«)(x- tj) 


1 

JV 


f>[x,  t=TjT+(k-l)T 


,j  =  l, 


(16) 


where  N  is  the  total  number  of  cycles  for  which  the  data  have  been 
sampled— a  minimum  of  200  for  the  authors’  data;  M  =  40  is  the  number 
of  samples  per  cycle;  and  T  is  the  period  of  the  applied  periodic  forcing. 
The  cyclic  modulation  of  the  turbulence  statistics  can  then  be  computed 
using  the  phase-locked  variance  operation: 

(q’2)(x,r.)=— ^-^{q[x,t  =  T.r+(fc-l)r]-(q)(x,r.)}2  ,j  =  1,  (17) 

jV  1  jc^1 


One-dimensional  profiles  of  phase-locked  statistics  were  then  computed  by 
streamwise  averaging  the  quantities  defined  in  Equations  16  and  17.  This 
procedure  is  valid  since  the  channel  flow  is  essentially  fully  developed  so 
that  gradients  in  the  streamwise  direction  are  expected  to  be  small. 

Instantaneous  flow  field  structure— two-point  spatial  correlations 

For  the  canonical  problem  of  fully  developed  channel/pipe  flow,  analytical 
solutions  exist  for  sinusoidal  time  variations  in  pressure  gradient  (Uchida 
1956;  Currie  1993)  or  volume  flow  rate  (Muntges  and  Majdalani  2002;  Ray 
et  al.  2005),  and  the  problem  is  well  understood.  In  the  turbulent  regime, 
experimental  data  that  describe  the  time-dependent  evolution  of  the 
phase-locked,  Reynolds-averaged  flowfield  statistics  have  been  provided 
(Hino  et  al.  1983;  Akhavan  et  al.  1991;  Jensen  et  al.  1989),  and  these 
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experimental  data  are  complimented  by  more  recent  reports  of  LES/DNS 
simulations  (Scotti  and  Piomelli  2001;  Costamagna  et  al.  2003;  Salon  et  al. 
2007).  Reynolds-averaged  information  is  a  useful  descriptor  of  turbulent 
fields  but  does  not  always  yield  the  needed  insight  into  the  fundamental 
dynamics  of  the  flow  that  is  given  by  the  identification  and  description  of 
spatially  and  temporally  coherent  structures,  which  have  been  extensively 
described  in  steady  turbulent  wall  flows  (Cantwell  1981;  Adrian  2007)  and 
have  been  found  to  contribute  significantly  to  the  production  and 
maintenance  of  wall  turbulence. 

The  description  of  the  time  evolution  of  coherent  structures  in  periodically 
unsteady  turbulent  wall  flows  (hereafter  referred  to  as  oscillating  flows')  is 
less  developed  than  for  the  steady  case,  but  a  picture  is  beginning  to  emerge 
based  largely  on  the  results  of  recent  numerical  simulations  (Scotti  and 
Piomelli  2001;  Costamagna  et  al.  2003;  Salon  et  al.  2007)  and  earlier  flow 
visualization  experiments  (Fishier  and  Brodkey  1991;  Sarpkaya  1993),  which 
must  be  verified  using  a  more  quantitative  experimental  description  based 
on  direct  and  spatially  correlated  measurement  of  the  velocity  field.  Earlier 
measurements  in  oscillating  turbulent  wall  flows  (Hino  et  al.  1983;  Akhavan 
et  al.  1991;  Jensen  et  al.  1989)  utilized  point  laser  doppler  velocimeter 
(LDV)  data,  from  which  structural  and  spatially  correlated  information 
cannot  be  extracted.  This  is  believed  to  be  the  first  PIV  measurements  in  a 
canonical  oscillating  wall  flow.  Earlier  research  provided  a  Reynolds- 
averaged  description  of  the  turbulent  field  along  with  a  brief,  qualitative 
presentation  of  the  instantaneous  2D  structure  of  the  flow.  A  portion  of  the 
PIV  data  set  is  analyzed  here  to  extract  quantitative  2D  structural  and 
spatial  correlation  information,  which  can  then  be  compared  to  the  present 
view  of  turbulence  structure  in  both  oscillating  and  steady  turbulent  wall 
flows. 

Statistical  imprints  of  the  instantaneous  structure  can  be  realized  using 
two-point  correlation  functions  of  the  velocity  fluctuations.  The  evolution 
of  two-point  self-correlation  functions  of  the  streamwise  (Run)  and  wall- 
normal  ( Rw )  velocity  fluctuation  fields  are  discussed  later  in  the  text.  Two- 
point  correlations  were  computed  from  phase-locked  velocity  fluctuation 
data  using 


Rij(r,y,y0,6) 


(ui{x,y0)uj(x  +  r,y)) 

°i{yo)°j{y) 


(18) 
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where  the  ij  subscripts  refer  to  the  streamwise  and  wall-normal 
directions,  the  Uk  represent  the  velocity  component  fluctuations  about 
their  respective  phase-locked  means,  the  ok  are  the  root  mean  square 
(rms)  fluctuations  of  the  uk,  and  the  (  }  brackets  represent  averaging  over 

the  x-coordinate  and  the  ensemble  of  200  phase-locked  PIV  realizations  at 
phase  angle  0.  The  correlations  were  computed  using  a  reference  height, 
y0  =  0.137//,  which  places  the  reference  point  in  the  log  layer  (when  the 
flow  is  turbulent)  for  the  majority  of  the  forcing  cycle. 


Instantaneous  flow  field  structure— conditional  velocity  fields 

The  two-point  spatial  correlation  fields  discussed  provide  a  description  of 
the  mean  size  and  orientation  of  the  turbulent  fluctuations  but  does  not  give 
further  specifics  on  the  details  of  the  structure  that  can  reveal  the  physical 
mechanisms  responsible  for  the  generation  of  correlated,  large-scale 
features.  In  steady  turbulent  wall  flows,  there  has  been  significant  recent 
evidence  to  link  boundary  layer  structure  and  the  production/maintenance 
of  Reynolds  shear  stress  to  the  existence  of  organized  groupings,  or  packets, 
of  hairpin  vortices  (Adrian  2007).  Several  questions  then  arise  regarding 
the  existence  and  history  of  hairpin  vortices  in  oscillating  wall  flows.  Are 
hairpins  present?  When  are  they  present?  How  do  hairpins  grow,  develop, 
and  dissipate  in  the  presence  of  strong  cyclic  acceleration,  deceleration,  and 
relaminarization?  To  begin  to  address  some  of  these  questions,  there  is  an 
extension  of  the  analysis  of  Christensen  and  Adrian  (Christensen  and 
Adrian  2001),  in  which  the  expected  value  of  the  velocity  fluctuation  field, 
given  the  presence  of  a  near-wall  swirling  strength  event— which  may  or 
may  not  arise  from  hairpin-vortex  motion— is  computed.  This  conditionally 
averaged  velocity  field  is  estimated  using  linear  stochastic  estimation,  which 
allows  a  computation  of  an  unbiased  estimate  of  the  average  based  on 
conditional  event  data  using  unconditioned  statistical  correlations  between 
velocity  and  swirling  strength.  The  proper  form  of  the  linear  estimate,  as 
given  by  Christensen  and  Adrian,  is 


Uj(x+r,y )  Ac;(x,y0) 


(19) 


From  Equation  19,  it  is  evident  that  the  expected  velocity  given  the 
condition  of  a  Ad  event  at  height  y0  from  the  wall  (note  that  the  i  in  Ad  is 
not  an  index  in  Equation  19)  can  be  computed  using  the  unconditional 
two-point  correlation  between  Ad  and  the  velocity  fluctuation  components. 
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Task  2-Create  data  base  of  shear-stress  time  histories  by 
demonstrating  and  applying  newly  developed  optical  sensing 
technology  within  SEAWOLF 

Identified  this  year  were  several  newly  developed  commercial  optical 
sensors  for  potential  rapid  shear-stress  measurement  applications  within 
the  SEAWOLF.  While  reliable,  PIV-derived  measurements  of  shear-stress 
time  histories  offer  limited  amounts  of  data.  Characterization  of  a  single 
wave-current  frequency  and  amplitude  operating  condition  generally 
requires  a  full  day  of  testing  and  an  additional  1-2  days  for  analysis  of  the 
PIV  data  for  the  extraction  of  wall  shear-stress  histories,  which,  to  date, 
has  created  a  limitation  of  shear-stress  results  for  only  eight  frequency/ 
amplitude  combinations.  For  the  purposes  of  obtaining  large  parametric 
data  sets  that  vary  frequency,  amplitude,  and  bed-roughness  conditions,  a 
more  direct  shear-stress  measurement  must  be  incorporated.  Floating- 
element  measurements  of  shear  stress,  based  on  load-cell  sensors,  provide 
a  direct  mechanical  force  measurement  that  can  be  calibrated  to  a  shear- 
stress  value.  These  mechanical  sensors  have  been  applied  for  testing  in 
SEAWOLF  with  mixed  results;  the  sensors  responded  to  time-varying 
shear  stresses  but  displayed  weak  signal  levels  that  were  easily  corrupted 
by  electrical  noise  and  suffered  from  time-dependent  drift  in  sensor 
voltage,  which  corrupted  the  measurement  as  well. 

As  an  alternative  to  mechanical  shear-stress  measurements,  there  was  an 
implementation  of  a  newly  developed  commercial  optical  sensor 
technology  that  (1)  offers  calibration-free  measurements,  (2)  improved 
signal-to-noise  characteristics,  and  (3)  has  been  proven  to  work  in  harsh 
and  dirty  environments.  Two  candidate  commercial  sensor  technologies 
from  Measurement  Science  Enterprise,  Inc.,  Pasadena,  CA,  were 
demonstrated.  They  were  the  microS  (Fourguette  et  al.  2004)  and  the 
miniLDV,  summarized  here. 

Micros,  a  micro-optical  sublayer  profilometer  and  shear-stress  sensor 

An  illustration  of  the  microS  sensor  platform  is  provided  in  Figure  9.  The 
sensor  contains  a  miniaturized  diode  laser  source  and  diffractive  optics 
that  project  a  coherent  fringe  pattern  into  the  near- wall  flow.  The  flow  is 
seeded  with  particles  (or  sediment  from  erosion  testing)  that  scatter  light 
from  the  fringe  pattern  that  is  sensed  by  a  fiber-coupled  detector  that 
images  a  measurement  volume  located  75-100  /nm  from  the  bed  surface, 
within  the  viscous  sublayer  where  the  velocity  profile  is  known  to  be  linear. 
The  probe  laser  and  detector  element  are  fully  integrated  into  the  compact 
sensor  housing  (Figure  9). 
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Figure  9.  Micros  optical  shear-stress  sensor. 


The  local  fringe  separation,  8,  was  designed  to  be  linear  with  the  distance 
from  the  sensor,  y,  given  by  5  =  k  y,  where  k  is  the  fringe  divergence  rate 
along  the  normal  to  the  sensor  surface.  As  seed  particles  flow  through  the 
linearly  diverging  fringes,  they  scatter  light  with  a  frequency,/,  which  is 
proportional  to  the  near-wall  velocity  and  inversely  proportional  to  the 
fringe  separation.  The  velocity  of  the  seed  particle  is  then  u  =  f8.  The 
frequency  simply  multiplied  by  the  fringe  divergence  yields 

—  =  fk  (20) 

y 

which  is  equal  to  the  wall  velocity  gradient,  provided  that  the  probe 
volume  is  located  within  the  viscous  sublayer.  The  bed  shear  stress  is  then 
calculated  from 


du  u  ri 

T  =  [X—=[i~=[lfk 

dy  y 


(21) 


where  (J  is  the  viscosity  of  water.  The  signal  conditioning  and  processing 
required  for  the  diverging  fringe  optical  sensor  is  identical  to  that  used  for 
a  standard  laser-Doppler  velocimetry  and  is  supplied  with  the  sensor  in  an 
integrated  package. 


This  commercial  sensor  can  be  easily  mounted  flush  with  any  bed  surface 
inside  an  enclosed  channel  facility  for  parametric  measurements  of  wall 
shear-stress  time  histories  with  rapid  turnaround.  It  is  easy  to  use, 
durable,  and  provides  reliable  data  in  the  presence  of  fouling  of  its 
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windows  by  seed  particles  and  sediment.  The  manufacturer  has  applied 
the  sensor  for  measurements  in  ocean  and  field  environments  where  it  was 
proven  to  be  robust. 

MiniLDV,  a  miniature  laser  doppler  velocimeter  (LDV) 

Most  LDVs  in  the  industry  are  table-top  systems  that  include  the  probe,  an 
argon-ion  laser,  fiber  manipulators  for  laser  alignment,  Bragg  cells,  and 
receiving  and  processing  enclosures.  Getting  results  quickly  requires  a 
user  with  experience  in  beam  alignment.  The  miniLDV  replaces  all  these 
components  with  a  probe  and  two  small  electronics  boxes,  decreasing  the 
sensor  system’s  overall  size  as  well  as  its  power  consumption.  Since  all  the 
sensitive  optics  are  enclosed  within  the  probe,  the  user  does  not  need  to 
spend  hours  setting  up  the  sensor,  leaving  more  time  to  perform  data 
collection. 

The  miniLDV  system  consists  of  a  transceiver  probe,  driver  electronics, 
bandpass  filter,  and  BP-miniLDV  burst  processor  acquisition  hardware 
and  software.  The  miniLDV  is  a  standalone  system  meaning  that  no  other 
lasers  or  optics  are  required.  For  use  it  is  only  necessary  to  place  the 
electronics  at  a  convenient  distance  from  the  acquisition  computer  and 
place  the  probe  at  the  desired  location.  The  burst  processor  software 
collects  data,  moves  the  probe  on  the  electronic  traverse,  and  presents  flow 
statistics.  The  process  of  measuring  a  profile  can  be  automated  with  the 
software  by  dictating  the  number  and  location  of  points  from  which  to 
collect  data. 

Conceptually,  an  LDV  system  consists  of  highly  coherent  light  split  into  two 
beams  and  sent  through  a  transmitter  that  crosses  the  beams  in  the  probe 
volume  to  create  a  fixed  interference  pattern  (Figure  10).  The  aberrations  of 
the  optics  must  be  low  so  that  the  fringes  of  the  pattern  are  straight 
throughout  the  probe  volume.  The  receiver  consists  of  light  collection  optics 
and  an  ultra-fast-response  photodiode.  When  a  particle  passes  through  the 
probe  volume,  it  passes  through  regions  of  destructive  interference 
(darkness)  and  constructive  interference  (brightness)  that  the  receiver  sees 
as  a  sinusoidal  wave  of  intensity.  Knowing  the  fringe  separation  (known 
from  calibration  of  the  given  transmitter),  the  frequency  of  the  sinusoid  is  a 
measure  of  the  speed  of  the  particle  through  the  probe  volume. 
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Figure  10.  Conceptual  operation  of  the  miniLDV  sensor. 


The  probe  volume  is  formed 
by  splitting  a  coherent  beam 
and  then  crossing  the  two 
halves  in  a  very  small  region 
to  generate  interference 
fringes. 


Probe  volume 


Row  direction 


As  the  particle  travels  through 
the  probe  volume,  it  reflects 
light  into  the  receiver  only 
when  it  crosses  a  region  of 
constructive  interference. 


Magnified  75  times 


The  frequency  of  the  sinusoid  is  measured 
(usually  by  FFT).  The  fringe  spacing  is 
known  (from  calibration),  thus  the  speed  of 
the  particle  can  be  calculated. 

speed  =  fringe  spacing  x/ 


As  seen  by  the  photodetector, 
the  light  signal  is  a  sinusoid 
(corresponding  to  the  fringes) 
modified  by  a  Gaussian 
envelope. 


) 


t 

The  sinusoid  is  isolated 
by  a  bandpass  frequency 
filter. 


The  miniLDV  includes  a  frequency  shifting  component  so  that  the  fringes 
in  the  probe  volume  are  moving  at  a  constant  velocity.  The  measurement 
then  becomes  one  of  the  relative  velocity  differences  between  the  particle 
in  space  and  the  fringes  in  the  probe  volume.  Thus  the  probe  can  specify 
flow  direction  as  well  as  measurement  speed. 

Task  3-Demonstrate  wave-current  bottom  shear-stress  time  histories 
representative  of  nature  in  SEAWOLF  by  manipulating  flume  and 
modulated  flow  conditions 

The  bottom  shear  stress  in  combined  current-wave  conditions  has  been 
extensively  studied  and  measured  (Grant  and  Madsen  1979;  Dibajnia  and 
Watanabe  1992;  Wikramanayake  and  Madsen  1994).  These  measurements 
have  been  used  to  develop  reliable  algorithms  for  bottom  shear-stress  time 
history  under  nonbreaking  waves.  The  various  algorithms  produce  fairly 
similar  results  and  in  general,  vary  as  to  the  amount  of  parameterization 
used  as  input.  For  example,  some  methods  require  near-bottom  velocity 
profile  time  series  as  input.  This  is  reasonable  for  laboratory  experiments 
but  is  not  appropriate  for  natural  systems  where  these  data  are  not 
available.  Available  data  generally  include  wave  height  and  period  as  well  as 
mean  current  velocity.  Several  studies  have  developed  algorithms  for  near¬ 
bottom  shear  stress  as  a  function  of  these  more  readily  available 
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parameters.  As  part  of  the  objective  to  replicate  naturally  occurring  wave- 
current  conditions  in  SEAWOLF,  these  parameterized  wave-current  shear- 
stress  methods  will  be  evaluated  against  calibrated  hydrodynamic  model 
output  and  used  to  program  piston  movement  and  head  difference  to 
produce  a  broad  range  of  user-specified  bottom  shear-stress  time  histories 
in  SEAWOLF. 

The  comprehensive  database  of  shear-stress  time  histories  collected  with 
the  optical  sensors  within  SEAWOLF  will  be  used  to  develop  a  means  to 
program  the  computer-driven  piston  movement.  A  detailed  parametric 
study,  which  includes  all  pertinent  combinations  of  wave  period, 
amplitude,  unidirectional  flow,  and  roughness  regimes,  is  necessary  to 
develop  an  understanding  of  how  each  parameter  affects  shear-stress  time 
history  in  SEAWOLF.  The  resolution  of  the  shear-stress  time  history 
database  (i.e.,  number  of  tests  performed)  will  determine  the  method  used 
for  filling  in  the  data  gaps  for  programming  the  SEAWOLF  flume 
operation.  High  resolution  in  the  database  lends  the  opportunity  for 
simple  data  fits  to  decipher  the  trends  associated  with  variance  in  the 
individual  parameters  affecting  shear-stress  time  history  while  a  low 
resolution  test  matrix  will  produce  a  more  uncertain  extrapolative  model. 
The  previous  near-wall  PIV  measurements  provided  high-quality  results 
that  revealed  the  nature  of  the  SEAWOLF  shear  stress  for  a  limited 
number  of  flow  conditions.  This  existing  database  will  be  significantly 
augmented  by  a  more  detailed  parametric  study  conducted  using  the 
miniLDV  sensor,  which  will  provide  shear-stress  data  in  a  much  more 
rapid  and  cost-effective  manner. 

Task  4-lnvestigate  uncertainty  in  shear-stress  measurement 
techniques 

The  wall  shear  stress  was  inferred  from  high- magnification  PIV 
measurements  of  the  velocity  profile  in  the  viscous  sublayer  using 
Equation  3.  The  impact  of  two  primary  causes  of  measurement  error  have 
been  considered:  (1)  uncertainty  in  the  viscosity,  //,  resulting  from 
temperature  variations  during  the  approximate  4-hour  duration  of  a  given 
measurement  run  and  (2)  the  confidence  interval  assigned  to  the  least- 
squares  estimate  of  the  slope  of  the  velocity  data,  du/dy,  in  the  viscous- 
sublayer  region. 


ERDC/CHL  TR-17-15 


29 


4  Results  and  Accomplishments 


Task  1-Determine  shear  stress  time  histories  and  turbulent  structure 
within  the  enclosed  SEAWOLF  channel  with  smooth  and  rough  walls 
under  modulated  flow 

Validation  of  PIV  shear-stress  measurements  in  a  steady  channel  flow 

The  validity  of  the  PIV  shear-stress  measurements,  obtained  using  the 
viscous-sublayer  approach  of  Equations  3  and  6b,  was  investigated  using 
data  from  steady  turbulent  flows  established  in  the  SEAWOLF  test  section. 
The  steady  turbulent  channel  flow  provides  an  excellent  test  bed  for  the 
validation  of  the  PIV  shear-stress  measurements  because  wall  shear  stress 
and  velocity  profiles  are  both  well  understood  in  steady  channel  flows. 
Wall  shear-stress  data  were  obtained  from  fits  to  both  the  linear  and 
logarithmic  laws  of  Equations  6b  and  10.  The  shear-stress  results  were 
additionally  compared  to  the  empirical  correlation  due  to  Blasius 
(Akhaven  et  al.  1991): 


-^%=0.0791Re71/4 
pU2  d 


valid  for  4000  <  Red  <  105 


(22) 


where  d  is  the  hydraulic  diameter  of  the  channel  and  U  is  the  mass-mean 
or  bulk  velocity,  U—l / 2 H  J  u{y)dy ,  with  H  representing  the  channel 

half  height  and  u(y)  the  PIV-measured  mean  velocity  profile. 

Sample  high-magnification  PIV  results  from  the  steady  channel  flow  at 
five  different  Reynolds  numbers  are  shown  in  Figure  6.  Mean  velocity 
profiles  were  computed  from  ensemble  and  streamwise  averaging  of  500 
statistically  independent  PIV  realizations,  as  shown  in  the  uppermost  plot 
in  Figure  8.  The  mean  velocity  data  were  normalized  using  inner-length 
scales  (Equation  6),  which  were  computed  using  shear-stress  values 
obtained  from  fits  of  data  in  the  logarithmic  layer  to  Equation  13.  At  all 
Reynolds  numbers  (based  on  bulk  velocity  and  hydraulic  diameter),  the 
velocity  profiles  display  good  collapse  to  a  universal  curve,  with  some 
scatter  observed  in  the  results  closer  to  the  wall,  which  is  indicative  of 
uncertainty  in  wall  location  and  small  differences  in  the  shear-stress 
values  computed  from  Equations  3  and  13. 
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The  resulting  wall  shear-stress  data  from  viscous-  and  log-layer  fits  are 
provided  as  a  function  of  Reynolds  number  in  the  lowermost  plot  in 
Figure  8.  The  experimental  results  are  plotted  alongside  empirical  values 
computed  from  Equation  (22).  Shear  estimates  from  linear  fits  to  the 
viscous  layer  are  within  -4%  to  -13%  of  both  Equation  (22)  and  results 
data  obtained  from  fits  to  the  log-layer  expression,  Equation  (13).  This  was 
seen  as  reasonable  level  of  accuracy  for  the  measurements  in  the 
oscillating  channel  flow.  Note  that  the  excellent  agreement  (within  a  few 
percent)  of  the  results  from  log-layer  fitting  and  Equation  (22)  is  likely 
fortuitous.  The  log-layer  fits  are  subject  to  uncertainty  in  the  exact  value  of 
the  constant  B  in  Equation  (13),  and  the  Blasius  correlation  has  additional 
uncertainty  that  is  unknown  to  the  authors  at  this  time.  There  is  no  reason 
to  believe  that  the  viscous  layer  fits  underpredict  the  true  value  of  the 
shear  stress,  and  it  is  the  agreement  of  all  three  methods  shown  in 
Figure  8  to  within  ~io%  that  is  encouraging. 

Phase-averaged  (see  below  for  a  discussion  of  phase-averaging 
procedures)  velocity  profiles  in  the  near-wall  region  for  test-matrix  Case  6 
are  presented  in  Figure  11.  The  profiles  are  presented  in  dimensional  form 
to  illustrate  the  high  spatial  resolution  in  the  near-wall  data  and  to  provide 
a  more  intuitive  link  to  the  slope  of  the  near-wall  profile,  shown  as  dark 
lines  in  Figure  12  and  the  wall  shear-stress  determined  from  Equation  (3). 
For  this  test  case,  51  mean  velocity  measurements  within  1  mm  of  the  test- 
section  wall  were  obtained,  from  which  at  least  5  and  as  many  as  11  points 
were  available  for  linear  fitting  within  the  viscous  sublayer. 

The  phase  modulation  of  the  wall  shear-stress  for  all  eight  cases  in  the  test 
matrix  of  Table  1  is  presented  in  Figures  12  (10  sec  wave  cases)  and  13  (5  sec 
waves).  For  each  case,  the  wall  shear  stress  (black  circles)  is  plotted  on  the 
left-hand  axis  and  the  volume  flow  rate  (white  circles)  on  the  right-hand 
axis.  The  shear-stress  results  were  obtained  from  the  slope  of  the  near-wall 
velocity  profiles,  similar  to  those  shown  in  Figure  11.  A  summary  of  the 
average  shear  stress  along  with  the  cycle  maximum,  minimum,  and  the 
amplitude  [(maximum  +  minimum)/2]  of  the  measured  shear-stress  waves 
is  provided  in  Table  2.  This  table  can  be  used  as  a  starting  point  for  the 
programming  of  the  SEAWOLF  facility  for  a  desired  shear-stress  history. 
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Table  1.  Summary  of  experimental  test  matrix.  Actual  metered  flow  rates  are  shown,  which  may  differ  slightly 

than  the  target  flow  rates  in  some  cases. 


Case# 

Max.  Flow 
(GPM)* 

Min.  Flow 
(GPM) 

Mean  Flow 
(GPM) 

Period  (seconds) 

Status 

1 

30 

-7 

11.5 

10 

Complete  wall  shear- 
stress  and  flowfield 
data.  Shaft  encoder 
used  as  a  time  base  for 

this  case.  Nonuniform  At 
between  data  points. 

2 

30 

-30 

0 

10 

Completed  high- 
magnification  PIV  and 
wall  shear  stress  for  20 
of  40  phases.  Low- 
magnification  images 
have  been  acquired  and 
are  currently  being 
analyzed 

3 

50 

-31 

10 

10 

Complete  wall  shear- 
stress  and  flowfield  data 

4 

50 

-50 

0 

10 

Complete  wall  shear- 
stress  and  flowfield  data 

5 

30 

-10 

10 

5 

Complete  wall  shear- 
stress  and  flowfield  data 

6 

30 

-30 

0 

5 

Complete  wall  shear- 
stress  and  flowfield  data 

7 

47 

-33 

7 

5 

Complete  wall  shear- 
stress  and  flowfield  data 

8 

50 

-50 

0 

5 

Complete  wall  shear- 
stress  and  flowfield  data 

*GPM  =  gallons  per  minute 


Table  2.  Summary  of  wall  shear-stress  results. 


CASE# 

1 

2 

3 

4 

5 

6 

7 

8 

Max.  flow  (GPM) 

30 

30 

50 

50 

30 

30 

47 

50 

Min.  flow  (GPM) 

-7 

-30 

-31 

-50 

-10 

-30 

-33 

-50 

Period  (SEC) 

10 

10 

10 

10 

5 

5 

5 

5 

Mean  shear  (Pa) 

0.65 

-0.04 

1.26 

-0.042 

0.53 

0.00 

0.81 

-0.15 

Max.  shear  (Pa) 

2.07 

2.54 

5.93 

6.13 

1.93 

2.24 

4.69 

7.72 

Min.  shear (Pa) 

-0.25 

-2.67 

2.96 

-6.50 

-0.39 

-1.94 

-2.38 

-7.51 

Amplitude  (Pa) 

1.16 

2.61 

4.45 

6.32 

1.16 

2.10 

3.53 

7.61 
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For  all  cases,  the  form  of  the  wall  shear-stress  history  is  significantly 
different  than  the  volume  flow  rate.  The  phase  lead  of  the  shear  stress 
relative  to  the  flow  is  not  particularly  significant,  which  is  evidence  of 
significant  turbulence  during  at  least  part  of  the  cycle.  As  a  check  on 
consistency  in  the  shear-stress  data,  note  that  the  cases  with  zero  mean 
flow  (2,  4,  6,  and  8)  display  shear-stress  waveforms  that  are  nearly 
symmetric  for  the  positive  and  negative  flow  portions  of  the  cycle. 
Measured  shear-stress  maxima  and  minima  for  the  zero-mean-flow  cases 
are  within  6%  for  three  of  the  four  zero-mean  cases,  with  a  maximum 
difference  of  15%  for  Case  6.  This  agreement  in  the  peak  values  of  the 
expected  symmetric  waveforms  is  an  initial  indicator  of  the  precision  of 
the  shear-stress  measurement. 

The  impact  of  the  cyclic  transition  between  laminar  and  turbulent  flow  is 
clearly  observed  in  the  shear-stress  results,  with  a  sharp  increase  in  shear 
stress  observed  with  transition  to  turbulence.  This  effect  is  especially 
apparent  in  the  data  for  Cases  1  and  5,  which  display  significant  shear- 
stress  jumps  during  the  later  stages  of  the  cycle  acceleration.  This  shear- 
stress  jump  was  well  documented  with  an  observed  flow  transition 
observed  in  the  PIV  results. 

For  each  plot  in  Figures  12  and  13,  an  additional  curve  labeled  “quasi¬ 
steady  (turbulent)”  has  been  added  to  the  graph.  This  curve  represents  the 
computed  value  of  the  wall  shear  stress  using  the  well-accepted,  steady 
turbulent  flow  result  from  Blasius,  given  earlier  in  Equation  (22),  which 
can  be  rearranged  to  solve  for  the  wall  shear  in  terms  of  fluid  properties 
and  the  measured  mass-averaged  velocity: 


r 


w 


0.0396pt77/4 


v1/4 


(23) 


where  vis  the  kinematic  viscosity  of  water.  The  cycle-resolved,  volume- 
averaged  velocity  can  be  evaluated  in  SI  units  from  the  measured  GPM 
data  using 


CPM 

U (m/s)=6.308  xlO  5  — -  (24) 

A 

where  A  is  the  cross-sectional  area  of  the  SEAWOLF  test  section.  The 
solid,  black  lines  in  Figures  12  and  13  were  computed  by  applying  the 
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steady  flow  correlation  in  Equation  (23)  but  using  the  unsteady  GPM  data 
to  obtain  the  quasi-steady  estimates  in  the  figures.  The  quasi-steady 
results  are  in  reasonable  agreement  in  their  predictions  of  average  shear 
stress  and  shear-stress  maximum  and  minimum  for  the  eight  cases  studied 
here,  so  this  formula  may  be  a  reasonable  way  to  interpolate  the  shear 
results  to  other  flow  conditions.  Equation  (23)  predicts  the  PIV-measured 
cycle-average  shear  stress  to  within  40%  and  the  maximum  or  minimum 
shear  to  within  27%,  with  some  cases  as  good  as  3%  agreement.  The 
correlation,  in  general,  overpredicts  the  PIV-measured  shear  stress. 

Cyclic  modulation  of  wall  shear  stress  and  turbulent  properties  in 
oscillating  channel  flow  within  the  SEAWOLF  flume 

A  test  matrix  of  eight  scenarios  (Table  1)  with  distinct  wave  current 
combinations  was  created  for  PIV  testing.  These  scenarios  adequately 
represent  the  range  of  operating  conditions  and  are  therefore  adequate  for 
CFD  calibration.  Low- magnification  PIV  studies  of  the  cycle-resolved  flow 
and  turbulence  in  the  SEAWOLF  channel  have  been  conducted  for  all 
eight  cases  in  the  test  matrix  of  Table  1.  The  consideration  of  just  one  of 
the  eight  cases  in  detail  should  be  illustrative  of  the  physics  of  flow  in  the 
flume.  Detailed  velocity  and  turbulence  data  for  test  matrix  Case  6  will  be 
discussed  in  this  report.  This  case  is  in  the  oscillating  flow  regime,  which 
is  characterized  by  zero  time-mean  velocity  (i.e.,  Umin  =  Umax).  Test  matrix 
cases  with  a  nonzero  time-mean  flow  are  characterized  as  pulsatile,  and 
the  presence  of  the  time-mean  velocity  component  leads  to  one  additional 
flowfield  parameter.  The  low-frequency,  large-amplitude  test  cases  in  the 
present  test  matrix  yield  pulsatile  flow  cases  which  appear  to  be  at  least 
qualitatively  similar  to  the  oscillating  test  cases  investigated.  Therefore, 
only  a  review  of  the  physics  of  the  oscillating  case  is  provided  here.  A  more 
detailed  review,  which  will  include  physics  relevant  to  pulsatile  flow,  will 
be  presented  in  upcoming  publications. 

The  reviews  by  Akhavan  et  al.  (1991)  and  Salon  et  al.  (2007)  identify  four 
regimes  for  wall-bounded  oscillating  Stokes  flows,  which  are 
parameterized  by  the  Stokes-thickness  Reynolds  number: 

ReS=U  o8/v  (25) 

where  U0  is  the  amplitude  of  the  cross-sectional  mean  velocity  oscillation 
and  <5  is  the  Stokes  thickness: 
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(26) 


vis  the  kinematic  viscosity  and  co  is  the  angular  frequency  of  the  velocity 
oscillation.  Specifically,  the  four  flow  regimes  are  (l)  laminar  flow  (L)  for 
Res  <  too;  (2)  disturbed  laminar  flow  (DL)  at  100  <  Res  <  500,  where  the 
velocity  modulation  displays  only  small  perturbations  about  the  laminar 
theory;  (3)  an  intermittently  turbulent  flow  (IT),  which  is  initiated  at  Res  ~ 
500-550;  and  (4)  a  full-cycle  turbulence  (T),  which  appears  for  higher  Res 
between  2600  and  3460,  based  on  the  inspection  of  the  shear-stress  time 
histories  in  Figure  9  of  Jensen  et  al.  (1989). 

All  the  PIV  experiments  were  conducted  in  the  IT  regime,  and  those 
discussed  here  had  Stokes-thickness  Reynolds  number  of  Res  =  1220, 
2033,  and  2875.  The  IT  regime  is  characterized  by  rich  physical  behavior, 
with  a  cyclic  relaminarization  as  a  result  of  very  strong  favorable  pressure 
gradient  during  the  cycle  acceleration  phase,  which  is  followed  by  a 
transition  back  to  turbulent  conditions  during  the  decelerating  or  more 
weakly  accelerating  portions  of  the  cycle.  As  Res  increases,  the  fraction  of 
the  cycle  for  which  the  flow  remains  turbulent  grows  until  full-cycle 
turbulence  is  eventually  achieved. 

A  small  number  of  quality  experiments  directed  at  turbulence  structure  in 
oscillating  wall-bounded  flows  have  been  reported  (Akhavan  et  al.  1991; 
Jensen  et  al.  1989;  Hino  et  al.  1983).  These  efforts  have  elucidated  the 
above-described  cyclic  turbulence-relaminarization-transition  process  and 
have  provided  quality  turbulence  data  utilizing  both  laser-Doppler  and 
thermal-anemometry  techniques  to  yield  pointwise  results,  but  which  do 
not  provide  any  instantaneous  picture  of  the  2D  structure  of  the  flow.  To 
date,  the  authors  are  unaware  of  any  reports  of  PIV  measurements  that  do 
provide  spatially  correlated  velocity  data,  and  the  present  measurements 
fill  this  gap.  The  PIV  results  described  here  have  also  yielded  extremely 
large  vector  ensembles,  arguably  resulting  in  better  converged  turbulence 
statistics  than  previously  reported.  The  PIV  implementation  reported  is 
additionally  unique  because  a  microscale  PIV  approach  to  a  macroscale 
problem  has  been  adopted  in  order  to  resolve  the  very-near-wall  structure 
and  because  the  dynamic  velocity  range  of  the  PIV  system  is  automatically 
adjusted  throughout  the  cycle  to  compensate  for  the  large,  periodic  swings 
in  velocity  magnitude. 
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Wall  shear-stress  behavior 

Phase-averaged  velocity  profiles  in  the  near- wall  region  for  Res=  1220  are 
shown  in  Figure  11.  The  profiles  are  presented  in  dimensional  form  to 
illustrate  the  high  spatial  resolution  in  the  near-wall  data  and  to 
demonstrate  the  quality  of  the  fits  to  velocity  data  in  the  viscous  sublayer, 
shown  as  dark  lines  in  the  figure.  For  this  test  case,  51  mean  velocity 
measurements  within  1  mm  of  the  test-section  wall  were  obtained,  from 
which  at  least  5  and  as  many  as  11  points  were  available  for  linear  fitting 
within  the  viscous  sublayer. 

Figure  11.  Representative  fits  to  near-wall  velocity 
profiles  recorded  in  an  oscillating  boundary-layer 
flow  at  Res=  1220.  Data  are  shown  in  dimensional 
units  and  plotted  on  linear  scales  to  illustrate  the 
effectiveness  of  near-wall  fitting. 
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Figure  12.  Phase-resolved  shear-stress  data  for  the  test-matrix  cases  with  10  sec  wave  period. 
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Figure  13.  Phase-resolved  shear-stress  data  for  the  5  sec  wave  period  cases. 
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The  phase-averaged  modulation  of  the  channel-averaged  velocity  is 


(u)= (u)(y)dy 


(27) 


wall  shear  stress,  turbulent  kinetic  energy,  and  production  rate  at 
Res=  1220,  2033,  and  2875  is  displayed  in  Figure  14.  The  channel- 
averaged  velocity  is  normalized  by  the  amplitude  of  the  average  velocity 
waveform,  Uo,  while  the  wall  shear-stress  results  are  presented  in  terms  of 
a  phase-averaged  friction  coefficient,  C/=  2<  Tw>/pU02.  The  time  histories 
of  the  in-plane  turbulent  kinetic  energy,  K,  have  been  averaged  across  the 
channel  height  and  normalized  according  to 


1 

2  HUl 


(u*)+(v'2)\dy 


(28) 


and  the  channel-averaged,  turbulence-production  rate  according  to 


2  HU30 


(29) 


where  the  Stokes  thickness,  5s,  has  been  selected  as  the  outer  length  scale. 


Figure  14.  Phase  modulation  of  the  bulk  velocity,  wall  shear  stress,  and  channel-averaged 

turbulence. 


ERDC/CHL  TR-17-15 


38 


For  all  cases,  the  form  of  the  wall  shear-stress  history  is  significantly 
different  than  the  nearly  sinusoidal  modulation  of  the  bulk  flow,  as  three 
terms  of  a  Fourier-series  representation  of  the  wall  shear-stress  data  are 
required  to  capture  the  main  features  of  the  shear-stress  oscillation.  At  all 
three  Reynolds  numbers,  the  first  harmonic  of  the  measured  shear-stress 
data  reveals  no  significant  phase  lead  of  the  wall  shear  relative  to  the  bulk 
velocity.  As  a  check  on  consistency  in  the  shear-stress  data,  note  that  the 
shear-stress  waveforms  are  nearly  symmetric  for  the  positive  and  negative 
flow  portions  of  the  cycle,  with  measured  shear-stress  maxima  and 
minima  within  15%. 

The  impact  of  a  cyclic  transition  between  laminar  and  turbulent  flow  is 
clearly  observed  in  the  modulation  of  <K>  and  <P>,  as  well  as  in  the  wall- 
friction  results.  At  each  Reynolds  number,  the  positive-velocity  half  cycle 
begins  with  <P>  at  essentially  zero  value  and  residual,  decaying  turbulent 
kinetic  energy  from  the  previous  half  cycle.  The  flow  is  relaminarized  in  the 
presence  of  a  strong,  favorable  pressure  gradient  during  the  acceleration 
phase,  with  an  explosive  transition  to  turbulence  observed  in  the  mid-to- 
late  acceleration,  dependent  on  the  value  of  Res.  The  transition  is  marked  by 
a  rise  in  <P>  from  near-zero  values,  and  this  rise  is  accompanied  by  a 
sudden  change  in  the  slope  of  the  friction-coefficient-vs.-time  curve.  This 
jump  is  particularly  apparent  at  Res=  2033  and  2875.  Note  that  production 
of  turbulence  ceases  for  some  fraction  of  the  cycle,  even  at  Res  as  high  as 
2875.  If  a  laminar  state  is  defined  as  one  where  <P>  is  at  near-zero  value, 
then  it  is  evident  that  the  laminar  portion  of  the  cycle  is  4096-50%  less 
duration  at  Res  =  2875  than  at  Res  =  1220,  based  on  inspection  of  the  phase- 
averaged  histories  of  <P>  in  Figure  14. 

For  comparison,  an  additional  curve  labeled  “quasi-steady,”  and  shown  in 
blue,  has  been  added  to  each  of  the  Cf  data  plots  in  Figure  14.  These  curves 
represent  the  computed  value  of  Cf  from  Equation  (14)  when  using  the 
phase-averaged  velocity  to  compute  the  hydraulic-diameter-based 
Reynolds  number,  Red.  The  quasi-steady  results  are  in  reasonable  agree¬ 
ment  in  their  predictions  of  average  shear  stress  and  shear-stress  extrema. 
Equation  (14)  predicts  the  PIV-measured  cycle-average  shear  stress  to 
within  40%,  and  the  maximum  or  minimum  shear  to  within  27%,  with 
some  cases  as  good  as  3%  agreement.  This  level  of  agreement  suggests  that 
the  quasi-steady  approach  may  provide  reasonable  engineering  accuracy. 
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Visualizing  the  bulk  flow  in  the  intermittently  turbulent  flow  (IT) 
regime  (introduction  to  flow  structure) 

The  case  at  Res=  1220  is  examined  here  in  more  detail  to  illustrate  the 
instantaneous  and  phase-averaged  structure  of  the  oscillating  channel  flow 
in  the  IT  regime.  The  cyclic  relaminarization  and  transition  process  is  well 
visualized  using  instantaneous  PIV  recordings  at  low  magnification,  as 
shown  in  Figures  15  and  16,  for  the  acceleration  and  deceleration  phases  of 
the  half  cycle,  respectively.  The  vector  fields  are  shown  across  the  full  depth 
of  the  channel,  and  the  instantaneous  value  of  the  globally  averaged 
streamwise  velocity  for  each  realization  has  been  subtracted  from  the  data 
to  better  reveal  the  structure  of  any  turbulence.  Vectors  are  displayed  using 
an  arbitrary  false-color  scale  as  an  indicator  of  the  velocity  magnitude;  the 
color  scale  indicates  the  difference  between  the  lab-referenced  and 
convective  velocities  and  is  red  for  slow-moving  fluid  and  blue  for  faster- 
moving  fluid  in  the  lab  frame  of  reference.  At  t/T  =  0.05,  the  velocity  field 
has  just  experienced  a  zero  crossing,  and  organized,  large-scale  structure  is 
observed  across  nearly  the  full  height  of  the  channel.  Production  of 
turbulence  is  negligible  at  this  time,  and  the  turbulent  kinetic  energy  is 
decaying  rapidly.  Relaminarization  then  abruptly  occurs,  with  the  vector 
fields  at  t/T  =  0.10  and  0.15  displaying  a  flat,  plug-flow  velocity  profile  in 
the  core  that  is  devoid  of  any  organized  large-scale  structure.  The  young 
laminar  boundary  layers  on  the  channel  walls  (observed  as  red-yellow  in 
Figure  15)  thicken  even  as  the  flow  is  accelerated.  By  t/T  =  0.25,  an 
instability  is  observed  near  the  channel  ceiling  and  in  the  core  flow; 
production  and  turbulent  kinetic  energy  begin  to  increase,  and  flow  reaches 
the  end  of  the  acceleration  phase  of  the  half  cycle. 

Instantaneous  vector  fields  for  the  decelerating  portion  of  the  half  cycle 
are  shown  in  Figure  16.  By  t/T  =  0.35,  when  the  bulk  flow  just  starts  to 
decelerate,  an  explosive  transition  to  turbulence  has  occurred;  large-scale 
structure  appears  in  the  core  of  the  channel,  and  turbulent  kinetic  energy 
and  production  are  at  their  cyclic  maxima.  Turbulence  persists  throughout 
this  period  of  strong  adverse  pressure  gradient,  through  t/T  =  0.40  and 
0.45,  while  the  rate  of  production  and  kinetic  energy  decline.  By  t/T  =  0.5, 
production  of  turbulence  has  ceased,  and  only  residual,  decaying,  but 
apparently  large-scale  structures  remain.  The  bulk  flow  then  changes 
direction  again,  and  this  process  repeats. 
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Figure  15.  Representative  instantaneous  velocity  vector  fields 
obtained  during  the  acceleration  phase  of  the  half  cycle  at 
Res=  1220. 


Figure  16.  Representative  instantaneous  velocity  vector  fields 
obtained  during  the  deceleration  phase  of  the  half  cycle  at 
Res=  1220. 
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Profiles  of  the  phase-averaged  streamwise  velocity  in  the  near-wall  region 
are  presented  in  terms  of  the  inner  variables  in  Figure  17.  The  profiles  are 
shown  for  Res=  1220,  and  these  data  have  been  scaled  using  the  phase- 
resolved  values  of  the  wall  shear  stress  shown  in  Figure  14  for  computation 
of  a  cycle-varying  friction  velocity  and  viscous  length  scale.  In  this  manner, 
a  comparison  to  the  law  of  the  wall  for  steady,  equilibrium  boundary  layers 
can  be  readily  made.  The  outstanding  spatial  resolution  of  the  high- 
magnification  PIV  results  readily  allows  an  observation  of  the  structure  of 
the  boundary  layer  from  the  logarithmic  region  to  well  inside  the  viscous 
sublayer.  In  each  profile,  the  viscous  sublayer  and  logarithmic  laws  (using 
k  =  0.41  and  B  =  5)  for  steady  boundary  layers  are  plotted  for  reference.  It 
is  evident  that  the  velocity  departs  significantly  from  the  logarithmic  law 
for  a  large  portion  of  the  cycle  so  that  fits  of  log-layer  data  cannot  be  used 
to  estimate  the  wall  shear  stress,  and  direct  probing  of  the  viscous  sublayer 
is  required.  The  half  cycle  begins  at  t/T  =  0.05  with  an  inflection  point  in 
the  profile,  and  the  velocity  data  lie  well  below  the  log  law.  By  t/T  =  0.1, 
the  flow  has  relaminarized,  and  logarithmic  behavior  is  not  observed  until 
t/T  =  0.3,  where  the  constant,  B,  in  a  log-law  fit  to  the  data  is  -6.5  for  k  = 
0.41.  This  logarithmic  behavior  persists  through  the  remainder  of  the  half 
cycle,  with  B  gradually  decaying  to  a  late-cycle  value  near  5.0  that  is 
typically  observed  in  steady  boundary  layers,  albeit  with  decreased  quality 
of  the  logarithmic  fit. 

Wall-normal  profiles  of  the  turbulent  stresses  at  Res=  1220  in  the  plane  of 
the  low-magnification  PIV  measurements  are  presented  in  Figure  18.  The 
low-magnification  PIV  results  reveal  the  cyclic  evolution  of  the  turbulence 
quantities  across  the  channel  half  height.  These  results  have  been 
normalized  using  the  amplitude  of  the  velocity  oscillation,  Uo,  and  the 
Stokes  thickness  as  velocity  and  length  scales,  so  that  the  data  are  more 
readily  compared  to  results  from  the  literature  (Jensen  et  al.  1989;  Salon 
et  al.  2007).  The  Reynolds  normal  stresses,  <u'2>  and  <v'2>,  have  been 
reported  in  terms  of  their  square  roots  (e.g.,  urms  =  ^(u'2) )  so  that 

variations  in  these  quantities  are  also  more  readily  evident.  The  Stokes 
thickness  is  further  indicated  by  the  dashed  lines  on  each  of  the  plots  in 
Figure  18. 
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Figure  17.  Near-wall  mean  velocity  profiles  at  Res=  1220 
displayed  in  inner  units.  The  viscous  scales  vary 
throughout  the  cycle  and  have  been  computed  using 
phase-averaged  values  of  the  wall  shear  stress. 


The  profile  of  the  streamwise  turbulence  intensity,  u'nns ,  is  flat  and  near 
zero  at  the  beginning  of  the  half  cycle.  From  there,  measurable  turbulent 
fluctuations  begin  very  near  the  wall  and  within  the  Stokes  layer  at  t/T  = 
o.i;  the  growth  of  the  streamwise  fluctuations  emerges  from  the  wall  and 
penetrates  to  the  core  of  the  channel  flow  by  t/T  =  0.3-0.35,  near  the  crest 
of  the  bulk  velocity  waveform,  after  which  the  strength  of  the  streamwise 
fluctuations  decay,  and  the  profiles  flatten  in  shape  until  a  flat  but 
measurable  streamwise  fluctuation  remains  just  after  the  zero-velocity 
crossing  near  t/T  =  0.5.  The  maximum  amplitude  of  the  streamwise 
fluctuations  remains  in  the  Stokes  layer  throughout  the  cycle,  with  its 
maximum  value  of  Urms/Uo  ~  0.17  at  t/T  =  0.25,  which  occurs  as  the  flow  in 
the  core  of  the  channel  is  just  beginning  to  transition  from  laminar  to 
turbulent  conditions  (Figures  15-16),  with  the  corresponding  jump  in  the 
wall  shear  stress,  as  seen  in  Figure  15.  The  streamwise  fluctuations  remain 
near  this  maximum  value  until  the  deceleration  phase  of  the  cycle  begins 
during  which  urms  decays. 
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Figure  18.  Phase-averaged  modulation  of  the  profiles,  in  outer  units,  of  the  turbulence  quantities  across  the 
channel  half  height  at  Res  =  1220,  as  acquired  using  low-magnification  PIV.  Streamwise  turbulence  intensity, 
Urms/Ua,  vertical  turbulence  intensity,  Vrms/Ua,  Reynolds  shear  stress,  <uV>/U<?\  and  rate  of  turbulence 
production,  &/Uo3  <uV>d<u>/dy.  Distance  from  the  wall  is  normalized  by  the  Stokes  thickness,  8&. 


The  behavior  of  the  streamwise  fluctuation  profiles  just  described  is 
similar  to  the  results  of  Salon  et  al.  (2007)  at  Res  =  1790  and  the  LDV 
measurements  of  Jensen  et  al.  (1989)  at  Res=  1000  and  1790.  Both  of 
these  earlier  studies  report  the  same  near-wall  generation  of  fluctuations 
and  subsequent  wave-like  penetration  of  these  fluctuations  into  the  core 
flow,  which  reaches  its  penetration  maximum  early  in  the  deceleration 
phase  at  t/T  ~  0.33-0.40,  as  compared  to  0.3-0.35  in  the  measurements 
at  Res  =  1220.  Jensen  et  al.  (1989)  report  maximum  streamwise 
fluctuations  of  Urms/Uo  ~  0.1  to  0.12  at  Res=  1000  and  1790  while  Salon  et 
al.  (2007)  report  maximum  values  of  Urms/Uo  near  0.15,  which  can  be 
compared  to  a  value  near  0.17  in  the  PIV  data. 

Vertical  velocity-fluctuation  profiles  in  Figure  18  reveal  that  peak  wall- 
normal  fluctuations  are  approximately  four  to  five  times  lower  than  the 
streamwise  fluctuations  observed  in  this  flow.  Peak  vertical  fluctuations 
are  in  the  vicinity  of  Vrms/Uo  ~  0.05,  compared  to  values  of  0.04-0.045 
measured  by  Jensen  et  al.  (1989)  at  Res  =  1000  and  1790,  and  computed 
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by  Salon  et  al.  (2007)  at  Res=  1790.  At  t/T  =  0.05,  the  measured  vertical 
fluctuations  display  a  monotonically  increasing  profile,  with  values  in  the 
core  that  are  comparable  to  the  streamwise  fluctuations  at  this  time.  The 
vertical-fluctuation  profiles  then  decay  in  the  channel  core  as  the  bulk  flow 
is  relaminarized,  with  a  second  inflection  point  observed  at  a  height 
corresponding  to  the  Stokes  thickness  discernable  at  t/T  =  0.10  and  0.15. 
The  vertical  fluctuations  then  proceed  outward  toward  the  channel  core  in 
the  same  wave-like  manner  displayed  by  Urms,  decaying  to  similarly  flat 
profiles  with  no  near-wall  maxima  by  the  end  of  the  forcing  cycle. 

Profiles  of  the  Reynolds  shear  stress,  and  rate  of  turbulence  production 
reveal  that  there  is  essentially  zero  production  early  in  the  half  cycle, 
production  begins  very  close  to  the  wall  at  t/T  ~  0.15,  and  the  turbulent 
shear-stress  and  production  emerge  from  the  near-wall  region  toward  the 
core  flow  in  the  same  wave-like  manner  discussed  above.  Production  of 
turbulence  is  confined  to  within  from  the  wall,  with  the  maximum  local 
value  of  production  occurring  within  the  Stokes  layer  at  t/T  =  0.25  and  the 
peak  integrated  production  at  t/T  =  0.35. 

Turbulence  structure  in  oscillating  channel  flow 

The  PIV  data  were  phase  locked  with  respect  to  the  piston  motion.  Data  at 
40  evenly  spaced  phase  angles  were  acquired,  with  200  velocity  realizations 
acquired  at  each  phase  angle  (8000  total  velocity  field  realizations).  For  the 
test  condition  to  be  discussed  here,  the  period  of  the  piston  forcing  was  5 
sec,  with  a  Uo  ±100  cm/sec  amplitude  velocity  wave  and  zero  time-mean 
flow  through  the  channel.  These  conditions  yield  a  Stokes  thickness  of  8 = 
1.23  mm  (from  Equation  26).  The  Reynolds  number,  from  Equation  19,  was 
2033,  which  is  in  the  fully  turbulent  regime  identified  by  Jensen  et  al. 
(1989).  Bulk-velocity,  wall  shear  stress,  and  turbulence  production/kinetic 
energy  waveforms  for  the  test  case  are  provided  in  Figure  19.  Under  these 
conditions,  the  shear  stress  varied  over  a  ±7.5  Pa  range,  such  that  the 
minimum  size  of  the  viscous  length  scale  was  11  /jm,  occurring  during 
periods  of  peak  wall  shear  stress.  The  spatial  resolution  of  the  PIV 
measurements  in  the  streamwise/wall-normal  plane  was  160  /an,  or  14.5 
minimum  viscous  wall  units.  This  represents  a  worst-case  spatial  resolution, 
as  the  PIV  vector  spacing  in  wall  units  correspondingly  improves  for  off- 
peak  shear-stress  conditions,  where  the  phase-resolved  value  of  the  viscous 
length  scale  is  larger. 
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Figure  19.  Phase-averaged  modulation  of  bulk  velocity  (top);  wall  skin  friction 
(middle);  in-plane  turbulent  kinetic  energy  and  turbulence  production  integrated 
across  the  channel  height  (bottom). 
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Instantaneous  flowfield  structure 

Representative  instantaneous  velocity  vector  realizations  are  shown  in 
Figures  20  and  21  for  the  accelerating  and  decelerating  portions  of  the  first 
half  cycle,  respectively.  The  vector  fields  are  plotted  for  the  lower  half  of  the 
channel  over  a  streamwise  distance,  x/H  of  1.8,  where  H  is  the  channel  half 
height.  A  streamwise  convective  velocity,  Uc,  of  85%  of  the  phase-locked 
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mean  at  the  channel  centerline  has  been  subtracted  from  the  vector  fields  to 
reveal  the  structure  of  turbulence  in  a  reference  frame  that  is  nominally 
moving  with  the  near-wall  flow.  Any  organized  vortical  structure  that  is 
traveling  at  Uc  will  appear  as  the  closed,  circular  vector  patterns  that  are 
associated  with  vortex  motion.  The  vectors  are  also  plotted  on  top  of  gray¬ 
scale  contours  of  swirling  strength,  Ad,  which  is  an  accepted  frame- 
independent  kinematic  identifier  of  vortex  motion  (Adrian  et  al.  2000a) 
that  marks  the  location  of  vortices  traveling  at  speeds  other  than  Uc. 

Instantaneous  turbulence  structure  during  the  acceleration  phase  (o°  to 
90°)  of  the  half  cycle  is  shown  in  Figure  20.  At  8°  and  170,  the  bulk  velocity 
is  near  zero,  and  reversals  of  the  bulk  flow  direction  and  sign  of  the  channel 
pressure  gradient  have  just  occurred.  At  this  time,  production  of  turbulence 
has  ceased  (Figure  15),  and  turbulence  from  the  decelerating  portion  of  the 
previous  half  cycle  persists  and  is  dominated  by  large-scale  vertical  up-  and 
downwash  motions.  Pockets  of  swirling  strength  are  observed  throughout 
the  channel.  By  350  and  530,  the  impact  of  strong  acceleration  is  felt 
throughout  the  flow.  Production  remains  essentially  zero  at  350,  and  weak 
swirling  strength  events  are  generally  limited  to  a  near- wall  region  for  y/H 
<  0.2;  the  flow  has  essentially  relaminarized.  At  530,  production  has  begun, 
with  some  vertical  ejection  of  swirling  motion  from  the  wall  observed  in 
some  fraction  of  the  realizations,  such  as  the  one  shown  in  Figure  20.  The 
acceleration  phase  ends  near  the  peak  of  the  velocity  wave,  where  the 
structure  is  typical  of  that  seen  at  710;  production  is  near  cycle-peak  values, 
an  organized  procession  of  strong  vortex  cores  reminiscent  of  hairpin  vortex 
signatures  (Adrian  et  al.  2000b)  appears  (more  evidence  below),  and  a 
logarithmic  phase-mean  velocity  profile  (Kearney  et  al.  2008)  is  established. 

Instantaneous  vector  and  swirl  fields  during  the  decelerating  portion  of  the 
half  cycle  are  displayed  in  Figure  21,  where  the  evolution  of  the  turbulence 
in  the  presence  of  a  strong  deceleration  is  visualized.  At  89°  and  107°,  the 
flow  is  in  the  same  quasi-equilibrium  state  as  at  710  above,  with  a  succession 
of  hairpin  vortex  cores  travelling  near  Uc  clearly  visualized  at  107°.  By  1250, 
a  weak  deceleration  has  initiated.  By  1610,  production  of  turbulence  is  still 
nonzero,  and  organized  vortical  structures  are  still  observed.  The  half  cycle 
ends  near  1790,  where  the  mean  velocity  profile  reveals  that  the  flow  in  the 
near-wall  region  has  already  changed  direction,  turbulence  production  has 
reverted  to  zero  levels,  and  very-large-scale  turbulence,  similar  to  that 
observed  near  the  start  of  the  half  cycle  at  8°,  reappears. 
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Figure  20.  Instantaneous  velocity  vector  measurements  during  the  acceleration 

portion  of  the  half  cycle. 


x/H 
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Figure  21.  Instantaneous  velocity  vector  measurements  during  the  deceleration  portion  of  the 

half  cycle. 
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Two-point  spatial  correlations 

The  time  evolution  of  the  two-point  correlation  fields  is  observed  in 
Figures  22  and  23,  where  contour  plots  of  Ruu  and  Rw  are  shown.  Note  the 
different  scales  on  the  horizontal  axes  in  Figures  22  and  23,  which  were 
used  because  the  streamwise  extent  of  the  Ruu  correlation  was  much  larger 
than  for  Rw.  For  quick  visualization  of  the  mean  size  and  orientation  of 
organized  turbulent  structure,  contours  for  R  values  higher  than  0.5  are 
filled  with  grayscale  colors,  with  lower-valued  contours  unfilled.  The  small 
ensemble  size  (200  PIV  realizations)  used  to  compute  the  correlations 
provides  only  marginal  convergence,  but  this  is  still  sufficient  to  visualize 
the  streamwise  evolution  of  these  quantities.  Larger  PIV  image  ensembles, 
while  desirable,  were  not  practical  as  a  result  of  the  unsteady,  phase-locked 
nature  of  the  experiments,  which  were  costly,  taking  6-8  hours  each  for 
setup  and  acquisition  of  200-250  forcing  cycles. 


Figure  22.  Two-point  spatial  correlation  of  the  streamwise  velocity  fluctuations  at  10  different  phase  angles  over 
the  first  half  cycle.  The  event  location,  yref,  is  located  at  y/H=  0.137. 
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Figure  23.  Two-point  spatial  correlation  of  the  wall-normal  velocity  fluctuations  at 
10  different  phase  angles  over  the  first  half  cycle.  The  event  location,  yref,  is 

located  at  y/H=  0.137. 
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Early  in  the  acceleration  phase  at  8°  and  170,  the  turbulent  fluctuations  are 
correlated  over  very  large  spatial  extents  with  the  Ruu  contours  exhibiting 
an  inclination  toward  the  upstream  direction  that  persists  from  the 
previous  half  cycle  where  the  bulk  flow  direction  was  reversed.  The  large 
spatial  extent  of  both  the  streamwise  and  wall-normal  velocity  correlation 
fields  is  associated  with  the  very  large  up-  and  downwash  structures 
observed  in  the  instantaneous  velocity  fields.  The  flow  is  relaminarized  at 
350,  where  the  extent  of  highly  correlated  Ruu  fills  the  channel  half  height 
with  diminished  size  observed  in  Rw,  which  can  be  interpreted  to  be  a 
result  of  cycle-to-cycle  jitter  in  the  observed  laminar  plug  flow  with 
essentially  zero  vertical  fluctuation.  By  530,  turbulence  production  has 
restarted,  and  a  highly  elongated  Ruu  is  observed,  with  no  discernable 
inclination  from  the  wall.  This  behavior  is  consistent  with  a  turbulent 
boundary  layer  in  the  near-wall  region  with  a  laminar  plug  flow  in  the 
channel  core.  The  elongated  Ruu  correlation  is  also  consistent  with  recent 
computational  observations  of  the  transition  process  in  oscillating  flows 
(Salon  et  al.  2007),  where  very  long,  streamwise-oriented  streaky 
structures  appear  in  the  near-wall  region  at  the  initial  onset  of  turbulence 
production.  From  710  to  1250,  the  Ruu  and  Rw  correlations  are 
quantitatively  similar  to  those  reported  by  Christensen  et  al.  (2004)  for 
steady  channel  flows  at  Rer=  1150.  In  the  710  to  1250  time  frame,  the  flow 
in  the  channel  is  near  the  peak  of  the  bulk  velocity  waveform  shown  in 
Figure  19,  with  minimal  temporal  change  in  pressure  gradient  and  a  near¬ 
constant  wall  shear  stress  resulting  in  a  value  of  ReT=  900.  The  Ruu  =  0.5 
threshold  is  ~o.gH  in  extent  with  an  upward  incline  in  Ruu  of  ~  io0-i2°  in 
the  data,  compared  to  12. 50  estimated  from  Christensens’s  (2004)  results. 
The  half  cycle  then  ends  with  a  growth  of  the  extent  of  both  Rm  and  Rw, 
and  the  next  half  cycle  then  begins  with  this  persistent  large-scale 
turbulence. 

Conditional  velocity  fields 

Linear  estimates  of  the  evolution  of  the  conditional  velocity  fluctuation 
fields  for  a  Ad  event  at  y0  =  0.137 H  (same  as  for  velocity  correlation  data 
above)  over  the  half  cycle  are  shown  in  Figure  24.  The  half  cycle  begins  with 
large-scale  residual  turbulence  that  persists  from  the  previous  deceleration 
phase;  no  appearance  of  the  in-plane,  hairpin-vortex  signature  of  the  type 
described  by  Adrian  et  al.  (2000b)  is  observed  until  the  late  acceleration 
phase,  with  the  first  evidence  of  hairpin  signatures  appearing  at  a  phase 
angle  of  62°  (not  shown  here),  where  acceleration  is  near  minimum  but 
production  of  turbulence  is  still  rising. 
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Figure  24.  Linear  stochastic  estimates  of  the  expected  velocity  field  given  a  Xa event  at  yref/H=  0.137 

and  x/H=  0. 


At  phase  angles  from  710  to  107°,  the  flow  near  the  peak  of  the  velocity 
waveform  has  achieved  an  essentially  quasi-steady  turbulent  state,  and  the 
conditional  velocity  fields  reveal  a  strong  vortical  circulation  at  the  swirl- 
event  location  (x  =  o ,y  =  y0)  with  other  features  indicative  of  the  presence 
of  hairpin  vortices  (Adrian  et  al.  2000b).  A  stagnation  line  in  the  reference 
frame  of  the  vortex  is  inclined  upward  in  the  direction  of  the  main  flow 


ERDC/CHL  TR-17-15 


53 


and  is  present  in  all  cases,  with  a  large  zone  of  uniform,  low-momentum 
fluid  underneath  the  stagnation  line.  This  shallow  inclination  angle  is 
similar  to  the  feature  observed  in  the  Ruu  contours  of  Figure  22  but  larger 
at  ~i6°.  The  vortex  cores  are  reminiscent  of  hairpin  vortex  heads  with  Q2 
events  upstream  of  the  vortex  and  associated  Q4  events  upstream  above 
the  stagnation  line  on  the  other  side  of  the  main  vortex.  Additional  weaker 
vortex  heads,  revealed  by  a  circular  vector  pattern  along  the  inclined 
stagnation  line,  are  observed  downstream  at  nominal  streamwise  spacings 
of  ~o.4 H,  suggesting  that  the  hairpins  are  indeed  traveling  in  packets 
convecting  at  approximately  the  same  speed. 

The  hairpin-vortex  signature  then  changes  during  the  deceleration  portion 
of  the  cycle.  From  126°  through  1610,  the  inclination  angle  of  the 
stagnation  line  is  reduced  although  the  presence  of  additional  vortices  at 
-0.4 H  spacing  persists,  until  the  end  of  the  half  cycle,  where  a  singular, 
isolated  vortex  is  observed  without  the  additional  features  of  hairpin 
vortex  packets  described  above. 

Task  2-Create  data  base  of  shear-stress  time  histories  by 
demonstrating  and  applying  newly  developed  optical  sensing 
technology  within  SEAWOLF 

Optical  shear-stress  sensors  demonstration 

A  representative  from  Measurement  Science  Enterprise,  Inc.,  performed  a 
demonstration  of  their  microS  and  miniLDV  sensors  in  January  2009. 
Since  these  measurements  were  taken  in  2009,  a  brief  summary  of 
findings  is  presented  here  while  more  detail  will  be  presented  in  upcoming 
reports. 

First,  the  microS  sensor  was  tested  under  unidirectional  flow  conditions. 
The  microS  reported  instantaneous  values  of  shear  stress  in  real  time  with 
the  software  provided.  This  allowed  for  analyzing  the  unidirectional  data 
for  the  average  shear  stress  as  well  as  turbulent  fluctuations.  Several 
different  unidirectional  flows  were  established  within  the  SEAWOLF,  and 
for  all  cases,  the  average  value  of  shear  stress  reported  by  the  microS  was 
within  10%  of  the  value  predicted  by  Blasius  (Equation  (14)).  This 
indicated  that  the  microS  worked  extremely  well  in  fully  developed 
unidirectional  turbulent  flow  cases. 
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Next,  the  microS  sensor  was  tested  under  oscillatory  flow  conditions 
within  the  SEAWOLF  flume.  Several  oscillatory  cases  were  tested,  but  in 
all  cases  the  quality  of  the  data  was  limited  by  low  data  rate  and  by  the 
limited  dynamic  range  of  the  shear-stress  sensor. 

The  miniLDV  sensor  was  found  to  be  a  much  better  suited  instrument  for 
the  measuring  of  wall  shear  stress  using  a  similar  single-point  estimate 
technique  than  the  microS  sensor  but  with  much  improved  data  rate  and 
sensor  dynamic  range.  The  miniLDV  sensor  is  mounted  externally  to  the 
SEAWOLF  test  section  and  is  insensitive  to  particle  fouling  on  sensor 
windows.  In  addition,  frequency  shifting  electronics  that  remove  the 
background  pedestal  from  the  Doppler  burst  signal  are  incorporated 
onboard  the  miniLDV,  which  provides  for  much  improved  sensor  dynamic 
range  that  is  a  must  for  the  wide  variations  in  shear  stress  observed  within 
the  oscillating  SEAWOLF  flow.  During  the  sensor  demonstration,  it  was 
experimentally  verified  that  the  miniLDV  velocity  signal  at  several  near¬ 
wall  locations  was  linearly  proportional  to  wall  shear  stress  for 
unidirectional  flow  in  SEAWOLF.  A  modified  miniLDV  sensor  from  MSE, 
with  ~30  ym  spatial  resolution,  was  purchased,  which  will  be  used  for 
detailed  parametric  studies  of  wall  shear  stress  in  Fiscal  Year  (FY)  09. 

Task  3-Demonstrate  wave-current  bottom  shear-stress  time  histories 
representative  of  nature  in  SEAWOLF  by  manipulating  flume  and 
modulated  flow  conditions 

Correlation  of  the  wall  shear-stress  results 

Akhavan  et  al.  (1991)  has  experimentally  shown  that  the  wall  shear-stress 
measurements  in  oscillating  pipe  flow  were  well  correlated  with  quasi¬ 
steady  applications  of  the  Blasius  formula  (Equation  (11))  during  turbulent 
portions  of  the  wave  cycle.  Akhhavan’s  laminar-phase,  shear-stress  data 
were  additionally  well  correlated  by  laminar  solution  for  oscillating 
channel  flow.  Comparisons  of  the  measured  wall  shear-stress  data  with 
quasi-steady  turbulent  predictions  were  presented  previously  in  this 
report.  Results  with  the  calculated  shear-stress  values  were  within  40%  of 
the  average  shear  and  27%  of  shear-stress  maximum  or  minimum. 

The  shear  stress  imparted  to  a  surface,  or  wall,  by  flowing  water  is  given  by 
Equation  (3).  With  LDV,  a  single-point  velocity  measurement  is  made  in  a 
region  near  the  wall  in  which  du/dy  is  well  approximated  by  a  linear 
gradient. 
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du 

dy 


\hjW  Uwall 
VlDV  —  ywall 


(30) 


When  the  no-slip  boundary  condition,  Uwaii  =  O,  is  applied,  and  the 
viscosity  to  be  a  constant  for  this  isothermal  water  flow  is  taken, 


t = constant  x- 


u 


LDV 


yLDV  Vwall 


(31) 


so  that  the  shear  stress  is  proportional  to  the  LDV-measured  velocity  when 
the  LDV  measurement  volume  is  at  a  fixed  distance,  ywv  -  ywaii,  from  the 
channel  wall. 

The  LDV  measurement  was  calibrated  for  shear  stress  by  performing 
measurements  in  a  steady  SEAWOLF  channel  flow,  where  the  wall  shear 
stress  was  previously  characterized  by  PIV  and  found  to  conform  to  well- 
known  results  for  turbulent  channel  flow.  This  results  in  an  approximately 
linear  calibration  curve,  as  shown  in  Figure  25,  where  the  measured  rvs. 
uldv  data  are  plotted  against  the  predicted  result  for  turbulent  channel 
flow  with  yLDV  -  ywaii  =  112  pm.  The  slope  of  the  measured  calibration 
curve,  as  shown  by  the  least-squares  fit  in  the  plot,  is  within  4%  of  the 
predicted  result.  Calibration  experiments  were  performed  before  each  LDV 
shear-stress  experiment,  and  linear  fits  of  the  type  shown  in  Figure  25 
were  used  to  convert  measured  LDV  velocities  to  wall  shear  stress.  Shear- 
stress  levels  encountered  in  unsteady  SEAWOLF  channel  flows  were 
observed  to  be  up  to  6  Pa,  or  twice  the  maximum  shear  present  during 
calibration.  The  linear  curves  were  extrapolated  to  these  higher  shear 
stresses  with  a  maximum  systematic  error  of  10%,  as  estimated  from  law- 
of-the-wall  velocity  profiles  for  turbulent  channel  flow. 
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Figure  25.  Measured  (symbols)  and  modeled  (solid  line)  LDV  shear-stress 
response  in  steady  channel  flow. 


Acquisition  of  initial  LDV  data  sets  and  comparison  of  LDV-  and  PIV- 
measured  shear  stresses 

Phase-averaged  shear  stresses  for  six  of  the  eight  SEAWOLF  operating 
conditions  monitored  by  PIV  have  been  measured  using  the  LDV  probe. 
The  results  for  5  sec  wave  periods  are  shown  in  Figure  26,  while  the  data 
for  10  sec  wave  periods  are  displayed  in  Figure  27.  Agreement  between  the 
PIV  and  LDV  data  is  reasonable,  with  shear-stress  amplitude  generally 
within  -30%  for  the  two  methods  and  similar  shear-stress  waveform 
shapes  measured.  The  LDV  probe  generally  yielded  higher  shear  stresses 
than  with  PIV.  In  some  cases,  the  measured  shear-stress  amplitudes  are 
within  10%.  Two  repeat  LDV  runs  are  plotted  in  Figure  28  for  the  case 
where  the  flow  rate  nominally  varies  between  +50  GPM  and  -50  GPM  at 
10  sec  wave  period.  The  repeatability  of  the  LDV  measurements  is 
excellent. 
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Figure  26.  Comparison  of  LDV-  and  PIV-measured  shear  stresses  for  5  sec  SEAWOLF  wave  period. 


Figure  27.  Comparison  of  LDV-  and  PIV-measured  shear  stresses  for  10  sec  SEAWOLF  wave  period. 
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Measurement  of  instantaneous  shear-stress  fluctuations  by  LDV 

One  major  advantage  of  LDV  probing  of  shear  stress  is  that  data  are 
acquired  on  an  instantaneous  basis,  versus  a  phase-averaged  basis  as  with 
PIV.  In  other  words,  only  the  average  shear-stress  wave  cycle  is  measured 
with  PIV  while  individual  wave  cycles  can  be  resolved  with  the  LDV 
sensor.  An  example  instantaneous  shear-stress  wave  form  is  shown  in 
Figure  28  for  a  case  with  maximum  and  minimum  flow  rates  of  30.2  and  - 
24.45  GPM  and  10  sec  wave  period.  Approximately  4  cycles  of  a  200+ 
cycle  data  set  are  shown  on  the  plot.  In  contrast  to  the  data  shown  in 
Figures  26  and  27,  where  only  the  phase-averaged  waveform  is  shown,  the 
instantaneous  shear-stress  wave  reveals  strong  turbulent  fluctuations 
superposed  with  the  phase  average.  Because  sediment  erosion  rate  is  often 
modeled  as  power-law  behavior, 


E=B  Tn 


(32) 


where  E  is  the  instantaneous  sediment  erosion  rate,  B  is  an  empirical 
constant.  The  exponent,  n,  is  an  empirically  determined  exponent  and  is 
frequently  found  to  be  greater  than  2  so  that  large  turbulent  fluctuations 
have  a  significant  effect,  with  an  impact  that  is  amplified  by  the  strong 
nonlinear  relationship  between  E  and  z.  LDV  measurements  of  the  kind 
shown  in  Figure  28  offer  the  possibility  of  quantifying  these  strong  erosion¬ 
enhancing  events,  while  PIV,  in  the  present  implementation,  cannot. 


Figure  28.  Instantaneous  shear-stress  wave  form  measured  using  LDV 

in  SEAWOLF. 
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Shear-stress  time  histories  of  the  type  shown  in  Figure  29  can  be  phase- 
decomposed  to  quantify  the  shear  stress  excursions.  Here,  the 
instantaneous  shear  stress  at  any  position,  x,  can  be  expressed  as  an 
average  waveform,  (r) ,  and  an  instantaneous  fluctuation,  r' ,  about  (r) : 

T(x,t)=(r)(x,0)+T’(x,t)  (33) 

Note  that  rand  r  are  functions  of  time,  while  (r)  is  a  function  of  the  cycle 

phase  angle,  0.  At  each  0  there  can  be  a  computation  of  an  rms  fluctuation 
of  r  about  (r) .  These  Trms  quantify  an  average  strength  of  the  shear-stress 

fluctuations  at  any  cycle  phase  angle,  and  the  shear-stress  fluctuations 
computed  from  the  full  200+  cycle  data  set,  of  which  a  subset  is  shown  in 
Figure  28,  is  shown  in  Figure  29.  The  baseline  shear-stress  fluctuation  is 
approximately  0.2  Pa.  At  these  cycle  times,  the  flow  is  relaminarized,  and 
the  baseline  fluctuations  are  likely  due  to  small  levels  of  jitter  in  the  cycle- 
to-cycle  flow  rate  waveform  and  by  measurement  error.  The  rms  shear 
stress  experiences  abrupt  increases  to  maximum  values  as  a  result  of 
transition  to  turbulence  and  decays  as  the  flow  is  decelerated.  The  flow 
changes  direction  at  a  phase-averaged  cycle  time  of  -5  sec  and  is 
relaminarized  again  before  the  process  is  repeated  in  the  negative  flow  rate 
portion  of  the  cycle. 

While  rms  shear-stress  does  a  good  job  of  quantifying  the  average  level  of 
turbulent  shear-stress  fluctuation  throughout  the  wave  cycle,  it  does  not 
capture  the  maximum  shear-stress  excursions.  To  address  this  issue, 
probability  densities  (histograms)  of  the  instantaneous  shear-stress 
measurements  have  been  computed.  A  representative  histogram  is  shown 
on  the  right-hand  side  of  Figure  30,  as  computed  from  the  data  set  of 
Figure  28.  Here,  maximum  shear-stress  excursions  of  5  Pa  (Figure  29, 
upper  right  corner)  are  observed;  it  is  the  tails  of  these  shear-stress 
histograms  that  can  lead  to  significant  unsteady  erosion  events. 
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Figure  29.  Representative  LDV  shear-stress  data  sets  for  two  different  flow  rate  conditions  at  10  sec  wave 

periods. 


Parametric  study  of  shear-stress  waveforms  generated  by  SEAWOLF 

The  LDV  sensor  in  production  mode  testing  has  been  applied  to 
investigate  the  shear-stress  waveforms  provided  by  SEAWOLF  under  a 
range  of  wave  and  current  flow  rates  and  for  wave  periods  of  5  and  10  sec. 
The  waveforms  were  generated  by  varying  the  following  controls: 

1.  The  programming  of  the  SEAWOLF  pistons  in  the  Anaheim  automation 
software;  four  different  programmed  piston  motions  were  run 

a.  “50_io”:  a  10  sec  wave  period  with  a  piston  stroke  of  8.17  in. 

b.  “50  5”:  a  5  sec  wave  period  with  a  piston  stroke  of  4.1  in. 

c.  “40_io”:  a  10  sec  wave  period  with  a  piston  stroke  of  6.54  in. 

d.  “40  5”:  a  5  sec  wave  period  with  a  piston  stroke  of  3.27  in. 

2 .  The  position  of  the  flow  metering  valve 

a.  Fully  open 

b.  450  to  the  horizontal 

c.  6o°  from  the  horizontal 
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3.  With  the  water  pump  that  provides  DC  or  steady  current  flow  on  or  off 

With  these  parameters  in  mind,  each  of  the  test  cases  was  assigned  a  unique 
name  so  that  it  was  obvious  what  the  SEAWOLF  control  settings  were  for 
each  experiment.  An  example  of  this  naming  convention  for  an  operating 
condition  running  the  “50  5”  SEAWOLF  piston  program  with  the  metering 
valve  position  at  450  to  the  horizontal  and  the  water  pump  off  is 

50_5_45_off 

All  data  sets  were  named  in  a  similar  manner,  with  a  total  of  22  different 
combinations  investigated  (Table  3). 


Table  3.  Summary  of  peak  phase-averaged  flow  and  shear-stress  data  in  SEAWOLF. 


Title 

T  (sec) 

Qmax  (GPM) 

Qmin  (GPM) 

t  max  (Pa) 

x  min  (Pa) 

50_5_45_on 

5 

50.8 

-29.6 

6.32 

-2.87 

50_5_60_on 

5 

52.46 

-41.4 

6.39 

-5.15 

50_5_open_off 

5 

31.3 

-29.37 

2.01 

-2.58 

50_5_45_off 

5 

43.4 

-40.84 

4.9 

-5.45 

50_5_60_off 

5 

48.7 

-44.97 

5.66 

-5.93 

50_10_open_on 

10 

57.15 

-17.44 

7.37 

-0.8323 

50_10_45_on 

10 

55.71 

-35.06 

6.98 

-4.46 

50_10_open_off 

10 

35.65 

-33.48 

3.71 

-4.09 

50_10_45_off 

10 

41.18 

-38.75 

4.72 

-5.15 

50_10_60_off 

10 

47.14 

-44.51 

5.36 

-5.78 

40_10_45_off 

10 

32.37 

-31.09 

3.38 

-3.77 

40_10_45_on 

10 

49.95 

-20.23 

5.57 

-1.23 

40_10_60_on 

10 

46.92 

-29.51 

5.2 

-2.96 

40_10_open_off 

10 

26.8 

-25.94 

2.26 

-2.37 

40_10_open_on 

10 

52.66 

-7.07 

5.53 

-0.1525 

40_5_open_on 

5 

48.26 

-3.2 

5.58 

0.0118 

40_5_45_on 

5 

44.71 

-27.04 

5.15 

-2.38 

40_5_60_on 

5 

43.25 

-33.91 

4.9 

-3.77 

40_5_open_off 

5 

23.42 

-22 

0.7448 

-0.8644 

40_5_45_off 

5 

29.33 

-28.09 

1.94 

-2.57 

40_5_60_off 

5 

37.51 

-35.13 

3.49 

-3.93 

50_5_open_on 

5 

52.79 

-12.39 

6.29 

-0.61 
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The  measurements  were  carried  out  with  the  LDV  probe  volume  positioned 
approximately  80  to  120  //m  from  the  floor  of  the  SEAWOLF  test  section.  At 
this  location,  the  probe  was  between  y+  of  4  to  10  viscous  wall  units  from  the 
wall  for  instantaneous  shear  stresses  ranging  from  1  to  6  Pa,  respectively. 
For  5  <y+  <  10,  the  probe  is  outside  the  viscous  sublayer  and  instead,  is  in 
the  lower  portions  of  the  buffer  layer,  where  small  departures  in  linearity  in 
the  velocity  profile  are  present.  Using  the  well-known  law-of-the-wall  for 
turbulent  velocity  profiles,  a  maximum  systematic  error  of  10%  for  a  phase- 
averaged  wall  shear  stress  of  6  Pa  is  estimated.  For  each  operating 
condition,  the  LDV  was  calibrated  for  its  shear-stress  sensitivity,  and  a 
curve  similar  to  the  one  shown  in  Figure  25  was  generated.  Oscillating  flow 
was  then  generated  in  SEAWOLF,  and  LDV  data  were  recorded 
continuously  for  a  period  between  20  and  30  minutes,  resulting  in  100  to 
300  wave  cycles  of  instantaneous  LDV  data.  The  calibration  was  repeated 
after  each  run  as  well  so  that  any  drift  in  the  velocity/shear-stress 
calibration  could  be  monitored  over  the  course  of  the  day’s  testing.  The 
calibration  was  then  applied  to  the  raw  LDV  data,  and  the  results  were 
phase  averaged  as  per  Equation  (5).  For  each  of  the  22  operating  conditions 
studied,  there  is  a  presentation  of  the  phase-mean  flow  rate,  phase-mean 
shear-stress,  phase-locked  rms  shear-stress  fluctuation  (the  fluctuation 
about  the  phase  mean),  and  a  10-wave-period  sample  of  the  instantaneous 
fluctuations. 

Sample  results  for  the  case  “40_io_open_off’  are  shown  in  Figure  30,  with 
the  data  for  the  remaining  21  operating  conditions  provided  in  Appendix  A. 
In  Figure  30,  the  phase-averaged  flowrate  (i.e.,  the  average  flow  rate  cycle) 
is  plotted  at  the  upper  left  with  the  corresponding  phase-averaged,  shear- 
stress  waveform  at  the  upper  right.  The  rms  value  of  the  phase-locked 
turbulence,  t  as  per  Equation  (5),  is  plotted  at  lower  left,  and  10  cycles  of 
the  instantaneous  LDV  data  are  at  lower  right  to  illustrate  the  degree  of 
turbulent  fluctuation  in  the  flow.  Note  that  hardware  for  locking  the  flow 
meter  and  LDV  probe  was  available  during  these  experiments,  and  the  zero 
time  of  the  flow  and  shear-stress  mean  and  rms  waveforms  was  adjusted  so 
that  the  positive-running  zero  crossing  was  approximately  at  time  =  o.  This 
was  done  to  provide  an  approximate  degree  of  correlation  between  the  flow 
and  shear-stress  data,  but  the  phase  relationship  between  these  two 
quantities  is  not  exact.  Transition  to  turbulence  with  the  associated 
profound  impact  on  shear  stress  can  be  identified  in  each  data  set  by  a 
sudden  spike  in  the  degree  of  scatter  in  the  instantaneous  shear-stress  data, 
which  is  accompanied  by  a  sudden  rise  in  the  phase-averaged  shear-stress 
and  rms  shear-stress  fluctuation. 
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Figure  30.  Phase-averaged  data  for  operating  condition  40_10_open_off. 


Results  for  all  22  cases  investigated  in  the  parametric  study  are  listed  in 
Table  3  in  terms  of  the  name  designation  of  the  data  set,  the  wave  period 
CT),  phase-averaged  minimum  and  maximum  flow  rates  (Qmin  and  Qmax ), 
and  phase-averaged  minimum  and  maximum  shear  stress  ( Tmin  and  Tmax). 
The  cycle  minimum  and  maximum  shear  stresses  were  plotted  against  the 
corresponding  minimum  and  maximum  flow  rates  for  all  22  data  sets, 
with  the  results  shown  in  Figure  31.  The  results  for  the  cycle  maximum 
(peak  or  positive  flow  rate)  are  plotted  in  blue,  with  the  data  for  the  cycle 
minimum  (trough  or  negative  flow  rate)  plotted  in  green  symbols.  Both 
sets  of  data  are  well  correlated  by  the  Blasius  shear-stress  equation  for 
steady  turbulent  channel  flow: 


r 


w 


0.0396 pU7/4 


(  \l/4 

V 


(34) 


and  this  level  of  agreement  is  consistent  with  the  quasi-steady  nature  of 
the  SEAWOLF  flow.  In  Equation  (34),  tw  is  the  wall  shear  stress,  p  is  the 
water  density,  U  is  the  volume-integrated  average  velocity,  d  is  the  channel 
hydraulic  diameter,  and  vis  the  kinematic  viscosity  of  water.  This  result 
shows  that  the  Blasius  equation  can  be  used  as  a  first-order  predictor  of 
the  amplitude  of  the  phase-locked  shear  stress  based  upon  the 
programmed  SEAWOLF  flow  rate  waveform. 
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Figure  31.  Cycle  minimum  and  maximum  phase-averaged  shear  stresses  correlated  with  their 
corresponding  minimum  and  maximum  phase-averaged  flow  rates. 


Task  4-lnvestigate  uncertainty  in  shear-stress  measurement 
techniques 

The  wall  shear  stress  was  inferred  from  high- magnification  PIV 
measurements  of  the  velocity  profile  in  the  viscous  sublayer  using 
Equation  (3).  The  impact  of  two  primary  causes  of  measurement  error  has 
been  considered:  (1)  uncertainty  in  the  viscosity,  //,  resulting  from 
temperature  variations  during  the  approximate  4-hour  duration  of  a  given 
measurement  run,  and  (2)  the  confidence  interval  assigned  to  the  least- 
squares  estimate  of  the  slope  of  the  velocity  data,  du/dy,  in  the  viscous- 
sublayer  region.  For  a  ±2  °C  variation  in  temperature  over  the  course  of  a 
given  data  run,  the  water  viscosity  was  found  to  vary  by  ±5%.  The  larger 
source  of  measurement  error  is  associated  with  the  confidence  interval 
associated  with  the  least-squares  estimate  of  du/dy.  In  general,  the  fitted 
slope  from  linear  fitting  of  the  sublayer  data  represents  the  true  value  of 

du/dy  within  an  interval  of  ±  t4c  ,  where  t  is  the  student’s  t  multiplier  for 
a  95%  level  of  confidence  and  C  are  elements  of  a  covariance  matrix  which 
measures  the  difference  between  the  measured  velocity  data  and  the 
velocities  predicted  from  the  least-squares  fit.  In  general,  the  confidence 
interval  improves  (get  smaller)  when  more  data  points  are  available  for 
fitting  and  as  the  signal-to-noise  in  the  PIV-measured  velocity  profile 
improves.  These  confidence  intervals  from  the  available  PIV  data  have 
been  calculated,  and  it  is  found  that  the  uncertainty  in  du/dy  is  generally 
within  ±20%  when  four  or  more  points  are  available  for  fitting  within  the 
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viscous  sublayer.  Based  on  these  results,  and  the  degree  of  scatter 
observed  in  the  measured  shear-stress  time  histories,  it  is  felt  that  the  data 
are  uncertain  to  within  ±20%  to  25%  or  better. 
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5  Concluding  Summary 

Task  1-Determine  shear-stress  time  histories  and  turbulent  structure 
within  the  enclosed  SEAWOLF  channel  with  smooth  and  rough  walls 
under  modulated  flow 

Cyclic  modulation  of  wall  shear-stress  and  turbulent  properties  in 
oscillating  channel  flow  within  the  SEAWOLF  flume 

A  hybrid-magnification  PIV  instrument  with  automated  dynamic-range 
adjustment  for  the  study  of  the  phase  modulation  of  the  wall  shear  stress 
and  flow  structure  in  an  oscillating  channel  flow  has  been  applied.  A  high- 
magnification  PIV  setting  was  utilized  to  obtain  superior  resolution  of  the 
near- wall  region,  with  an  in-plane  spatial  resolution  ranging  from  0.2  to  1.6 
viscous  wall  units,  which  enabled  high-resolution  observation  of  the  phase 
modulation  of  the  viscous  and  logarithmic  regions  of  the  velocity  profile  in 
oscillating  flow.  At  Res  =  1220,  it  was  observed  that  relaminarization 
inhibits  the  formation  of  a  logarithmic  velocity  profile  until  well  into  the 
later  portions  of  the  velocity  half  cycle.  Fitting  of  late-cycle  velocity  data  at 
Res  =  1220  to  a  logarithmic  curve  with  k=  0.41  reveals  that  the  constant  B 
decays  from  -6.5  at  the  inception  of  logarithmic  behavior  to  a  value  near 
5.0  at  the  end  of  the  half  cycle.  Fits  to  the  mean  velocity  data  in  the  viscous 
sublayer  were  used  to  obtain  reliable  wall  shear-stress  measurements, 
which  were  within  ~10%-15%  of  shear-stress  estimates  obtained  from  an 
accepted  empirical  correlation  and  fitting  to  log-layer  velocity  data  in  a 
steady  flow  test  case.  A  low-magnification  PIV  setting  was  used  to  reveal  the 
flow  structure  across  the  full  channel  height,  revealing  the  instantaneous 
flow  structure  during  relaminarization  and  transition  at  Res=  1220  and 
quantifying  the  phase-evolution  of  the  channel-integrated  turbulent  kinetic 
energy  and  production  rate  at  Res  =  1220,  2033,  and  2875.  The  results 
quantify  the  impact  of  transition  and  relaminarization  on  wall  shear  stress 
and  quantify  the  decrease  in  the  duration  of  the  laminar  phase  of  the  cycle 
by  about  a  factor  of  2  for  Res  from  1220  to  2875. 

Turbulence  structure  in  oscillating  channel  flow 

A  detailed  analysis  of  the  structure  of  turbulence  in  an  oscillating  channel 
flow  at  Res=  2033  using  PIV  data  in  the  streamwise/ wall-normal  plane  has 
been  presented.  To  the  authors’  knowledge,  this  is  the  first  structural 
information  for  an  oscillating  flow  derived  from  direct  velocity  measure¬ 
ments.  The  phase  evolution  of  the  two-point  spatial  correlation  functions, 
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Ruu  and  Rw,  along  with  linear  stochastic  estimates  of  the  velocity  field  in  the 
presence  of  a  swirling  strength  event  in  the  log  layer  are  presented.  The 
results  reveal  a  flowfield  that  is  dominated  by  very  large  residual  turbulence 
events  in  the  early  phases  of  the  cycle  acceleration  phase,  which  is  followed 
by  a  brief  period  of  relaminarized  flow  in  the  presence  of  a  strong 
acceleration.  Turbulence  re-emerges  in  the  later  portions  of  the  cycle 
acceleration  phase,  where  the  weakly  accelerated/decelerated  flow  in  the 
vicinity  of  the  peak  bulk  velocity  attains  a  quasi-steady  turbulent  state  with 
evidence  of  hairpin-like  structures  resembling  those  reported  in  steady 
turbulent  wall  flows.  These  hairpin-like  vortices  are  observed  to  persist 
through  much  of  the  deceleration  phase  of  the  half  cycle,  vanishing  only  just 
before  the  end  of  the  half  cycle  where  flow  reversal  and  a  change  in  the  sign 
of  the  acceleration  occur  in  the  near-wall  region. 

Task  2-Create  data  base  of  shear-stress  time  histories  by 
demonstrating  and  applying  newly  developed  optical  sensing 
technology  within  SEAWOLF 

Optical  shear  stress  sensor  demonstration 

As  an  alternative  to  laborious  PIV  shear-stress  measurements,  two  newly 
developed  commercial  optical  sensor  technologies  were  tested,  the  microS 
and  miniLDV,  which  (l)  offer  calibration-free  measurements,  (2)  rapid 
interrogation  of  near  bed  velocities  for  wall  shear  stress,  and  (3)  have  been 
proven  to  work  in  harsh  and  dirty  environments.  The  miniLDV  significantly 
outperformed  the  microS  in  initial  testing  and  therefore  was  implemented 
within  the  SEAWOLF  facility. 

Task  3-Demonstrate  wave-current  bottom  shear-stress  time  histories 
representative  of  nature  in  SEAWOLF  by  manipulating  flume  and 
modulated  flow  conditions 

Comparisons  of  the  measured  wall  shear-stress  data  with  quasi-steady 
turbulent  predictions  were  presented.  Results  with  the  calculated  shear- 
stress  values  were  within  40%  of  the  average  shear  and  27%  of  shear-stress 
maximum  or  minimum.  While  this  may  seem  high,  shear-stress  fluctuations 
in  natural  systems  are  typically  of  a  similar  range.  A  total  of  22  parametric 
flow-rate  period  and  amplitude  experiments  were  conducted  to  generate  an 
empirical  data  base  that  can  be  used  to  program  the  SEAWOLF  facility  for  a 
desired  shear-stress  waveform  to  be  used  in  erosion  testing  of  sediment  core 
samples.  The  parametric  data  were  generated  for  wave  periods  of  5  and  10 
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sec  with  flow-rate  amplitudes  ranging  from  5  to  60  GPM  for  cases  with  and 
without  unidirectional  current  flow.  The  results  reveal  the  complex 
periodically  unsteady  laminar-to-turbulent  transition  followed  by 
relaminarization  processes  in  SEAWOLF,  with  peak  shear-stress  amplitudes 
primarily  occurring  in  the  turbulent  regime  and  well  correlated  to  peak  flow 
rates  using  the  well-know  Blasius  correlation  for  steady  turbulent  channel 
flow.  Results  indicate  that  the  Blasius  equation  (Equations  (19)  and  (20)) 
can  be  used  as  a  first-order  predictor  of  the  amplitude  of  the  phase-locked 
shear  stress  based  upon  the  programmed  SEAWOLF  flow  rate  waveform. 
This  data  set  can  now  be  used  in  a  practical  fashion  to  provide  shear-stress 
time  histories  for  erosion  testing. 

Task  4-lnvestigate  uncertainty  in  shear-stress  measurement 
techniques 

Consideration  has  been  given  to  the  impact  of  two  primary  causes  of 
measurement  error:  (1)  uncertainty  in  the  viscosity,  //,  resulting  from 
temperature  variations  during  the  approximate  4-hour  duration  of  a  given 
measurement  run  and  (2)  the  confidence  interval  assigned  to  the  least- 
squares  estimate  of  the  slope  of  the  velocity  data,  du/dy  in  the  viscous- 
sublayer  region.  For  a  ±2  °C  variation  in  temperature  over  the  course  of  a 
given  data  run,  the  water  viscosity  was  found  to  vary  by  ±5%.  The  larger 
source  of  measurement  error  is  associated  with  the  confidence  interval 
associated  with  the  least-squares  estimate  of  du/dy.  Calculation  of 
confidence  intervals  from  the  available  PIV  data  demonstrate  that  the 
uncertainty  in  du/dy  is  generally  within  ±20%  when  four  or  more  points  are 
available  for  fitting  within  the  viscous  sublayer.  Based  on  these  results,  and 
the  degree  of  scatter  observed  in  the  measured  shear-stress  time  histories,  it 
is  felt  that  the  data  are  uncertain  to  within  ±20%  to  25%  or  better. 
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Appendix  A:  Compilation  of  Shear  Stress  Data 


Figure  Al.  Shear-stress  results  for  operating  condition  50_5_45_on. 


Figure  A2.  Shear-stress  results  for  operating  condition  50_5_60_on. 
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Figure  A3.  Shear-stress  results  for  operating  condition  50_5_open_off. 
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Figure  A4.  Shear-stress  results  for  operating  condition  50_5_45_off. 
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Figure  A5.  Shear-stress  results  for  operating  condition  50_5_60_off. 


Figure  A6.  Shear-stress  results  for  operating  condition  50_10_open_on. 
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Figure  A7.  Shear-stress  results  for  operating  condition  50_10_45_on. 


Figure  A8.  Shear-stress  results  for  operating  condition  50_10_open_off. 
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Figure  A9.  Shear-stress  results  for  operating  condition  50_10_45_off. 


Figure  A10.  Shear-stress  results  for  operating  condition  50_10_60_off. 
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Figure  All.  Shear-stress  results  for  operating  condition  40_10_45_off. 


Figure  A12.  Shear-stress  results  for  operating  condition  40_10_45_on. 
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Figure  A13.  Shear-stress  results  for  operating  condition  40_10_60_on. 


Figure  A14.  Shear-stress  results  for  operating  condition  40_10_open_off. 
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Figure  A15.  Shear-stress  results  for  operating  condition  40_10_open_on. 


Figure  A16.  Shear-stress  results  for  operating  condition  40_5_open_on. 
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Figure  A17.  Shear-stress  results  for  operating  condition  40_5_45_on. 


Figure  A18.  Shear-stress  results  for  operating  condition  40_5_60_on. 
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Figure  A19.  Shear-stress  results  for  operating  condition  40_5_open_off. 
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Figure  A20.  Shear-stress  results  for  operating  condition  40_5_45_off. 
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Figure  A21.  Shear-stress  results  for  operating  condition  40_5_60_off. 


Figure  A22.  Shear-stress  results  for  operating  condition  50_5_open_on. 
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